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I.  I  INTRODUCTION 

I 

jf 

During  tlie  past  year  the  Liiser  Kinetics  Branch  has  been  studying  the  B^II  -  X^Il 
transition  in  nitrogen-sulfide.  We  liave  been  particularly  interested  in  this  transition 

4'.  -  t*  . 

because  we  think  the  B  states  1/t  second  lifetime  and  large  h'ranck-Condon  shift  make  it  a 
good  short  wavelength  chemical  laser  candidate.  As  part  of  our  study  we  have  determined 
till'  rotational  and  vibrational  populations  of  the  excited  molecule  produced  in  an  S+N3 
NS(B)  +  No  reaction.  To  do  this  accurately  it  was  necessary  to  write  a  spectrum  code, 
d'liis  document  outlines  the  theory,  explains  how  to  run,  and  provides  a.  listing  of  the 
spectrum  code  and  accompanying  files.  It  should  bo  noted  that  even  though  this  code  was 
written  siiecifically  for  the  NS(n''^ll)  -  NS(X2J1)  transition,  it  is  easily  modified  to  do  a 
variety  of  other  molecules  and  transitions.  j  4.  " 

II.  NS  SPECTRUM  CODE  TIIEOHY 


The  code  is  based  on  the  quantum  mechanics  of  Schrbdinger.  In  the  next  few  pages, 
I'll  outline  the  most  significant  relations  loading  up  to  the  oloctronic-ro-vibrational 
intensity  relation  which  gives  the  intensity  of  a  spectrum  line  as  a  function  of  wavelength. 


The  perturbed  Schrddinger  eciuation  is  given  by. 


(2.1) 


where  nik  is  the  mass  of  the  k*'*'  particle,  \  is  the  wave  function;  Xk,  Yk,  Zk  are  the 
coordinates  of  the  k'*’  particle,  E  is  the  total  energy,  V  is  the  i)oteiitial  energy,  and  a;,\  is 
the  i)orturbing  term.  'I'he  interaction  of  an  electromaguetic  wave  (with  an  electric  vector 
E)  with  an  atomic  system  can  be  a.|)proxi mated  as  the  interaction  of  E  with  the  dipole 
moment  (M)  'of  the  atomic  system.  'I'liis  interaction  enters  into  Icquation  2.1  as  u)  =  E-M, 
where  M  has  components  Mx=^'kXk,  My=SekYk,  M^-iijkZk,  and  Ok  is  the  change  on  the 
kih  particle. 

In  solving  the  Schrbdinger  eciuation  for  an  electromagnetically  disturbed  atomic 
system,  Pauling  and  Wilson'  found  that  if  the  system  was  originally  in  a  state  of  energy. 


E,„  there  is  a  non-'/ero  probability  of  niidiiig  it  in  a  state  (^f  energy,  K,,,,  il'  the  emission  of  a 
photon  of  energy  hc/^nm  =  I'mi  ~  Em  allowed.  Specifically,  the  i)robability  of  a  transition 
between  two  states  n  and  ni  is  proportional  to  the  square  of  the  matrix  elements  of  tlic 
electric  dipole  moment  (R""')- 

The  comiwnents  of  are  defined  as  follows: 


Rx""’  =  J  Xn*MxX'tn‘l7’, 

Rynm=  J  x,,*MyXm(}r, 


(2.2) 

(2.3) 


R^nm  =  J.\n*M2,Vmdr,  (2.4) 

where  Xn,  Xm  are  the  eigenfunctions  of  the  two  states,  and  the  integral  is  over  the  total 
configuration  space. ^ 

The  intensity  of  a  spectral  line  for  a  transition  from  state  Xn  to  Xm  is  given  by: 


leu/’"'  ~  Nnhc/>'iiiiiAniii, 


(2.r>) 


where  N„  is  the  number  of  atoms  or  molecules  in  the  initial  stale,  iicz/nm  =  En-Em  =  energy 
emitted  in  transition,  A„m  is  the  ICinstein  transition  probability  of  spontaneous  emission 
and  according  to  wave  mechanics  is  related  to  the  matrix  elements  of  the  electric  dipole 
moment  (as  was  previously  shown  for  stimulated  emission)  as  follows: 

Ann  =  for  non-degenerate  levels. 

Anm  =  ^  for  levels  that  are  d,>  fold  degenerate.  (2.G) 


Combining  Equations  (2.5)  and  (2.6)  gives  the  intensity  of  the  .‘spectrum  line  in 
emission  as: 

However,  our  experimental  detectors  only  count  the  number  of  i)hotons  in  a  given 
frequency  range.  The  number  of  counts  at  a  given  frecpiency  can  be  obtained  by  dividing 


2 


Equation  (2.7)  by  liie  energy  associatetl  with  each  photon  at  that  frequency.  Doing  this 
yields  the  photon  count  as  a  function  of  freriuency; 


T  nm  _  M  ^ I  i '"k  !  ^  lno\ 

~W~  11,7 

The  number  of  molecules,  N,„  in  the  initial  state  depetid  on  the  electronic, 


vibrational,  and  rotational  levels  available,  as  well  as  the  energy  required  to  populate  each 
level.  In  the  following,  we  will  determine  the  energy  required  to  populate  a  given 
clcctronic-ro-vibrational  state  as  well  as  tlic  levels  available.  For  equilibrium  systems  we 
know  these  states  to  be  populated  according  to  a  Boltzmann  distribution.  This  gives  the 
number  of  molecules  in  each  state  once  we  know  the  total  energy  available  to  the  system. 
11.1.  DETFBMl.NA'nON  OF  ENEIUIY  LEVELS 


We  know  that  the  energy  of  a  molecule  is  due  to  lh('  motion  of  tlu*  ('l('ctrous  about 
the  nuclei  (electronic  eiungy),  vibration  of  tiie  nuclei  (vibrational  energy),  and  rotation  of 
the  nuclei  about  the  center  of  mass  (rotational  energy).  We  first  will  consider  the 
rotational  energy.  If  we  api)roximate  the  molecule  as  a  rigid  rotator  and  substitute 
m  =  ^  =  and  V  =  0  into  the  Shrbdinger  eriuation  (E(|uation  2.1),  we  obtain  a 

solution  for  the  rotational  energy  levels; 


where  Mj,  M2  are  the  nuclear  masses,  and  1  is  the  moment  of  inertia  of  the  rotator;  J  refers 


to  the  rotational  level  in  question. 

Now  in  spectroscopy  it  is  customary  and  convenient  to  write  all  energies  in  terms  of 
wave  numbers.  The  wave  number  is  obtained  by  dividing  the  energy  by  he.  So,  in  wave 


numters  E(|uation  (2.9)  becomes 
E  I  1.  I  /  I  I  1  \ 


U  =  (enr')  (2.10) 

B  =  I  is  called  the  rotational  constant.  'I'lic  wave  mimber  {or  energy)  associated 


with  a  transition  from  one  rotational  level  to  another  is: 
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i/=  Ff”)  -F(J"), 


(2.11) 


where  .1'  always  refers  to  the  upper  rotational  level  and  J"  refers  to  the  lower. 

It  is  well  known  that  molecules  not  only  rotate  alx)ut  their  center  of  mass,  the  nuclei 

vibrate  with  respect  to  one  another.  For  a  diatomic  molecule,  the  molecule  vibrates  along 

an  imaginary  line  joining  them  and  passing  through  their  respective  centers.  The 

coulombic  force  holding  the  nuclei  together  goes  as  -kx.  The  potential  energy  associated 
kx'^ 

with  this  is  V  =  —:y-  or  that  of  a  harmotiic  oscillator.  If  one  substitute.s  this  into  the 
Schrodinger  equation  and  solves  for  the  energies  of  the  vibrational  levc'Is,  the  following 
wave  numbers  (or  energies)  are  ol)taine(l: 

^  =  (l(v)  =  a;(v  +  i),  (2.12) 

where  u  refers  to  the  fundamental  fre(|uency  of  vibration  lx>tween  the  nuclei,  and  v 
corrcsfX)nds  to  the  allowed  vil)rational  levels.  The  wave  number  associated  with  a 
transition  from  one  vibrational  level  to  another  is: 

t/=  G(v')-G(v"),  (2.13) 

where  v'  refers  to  the  upper  vibrational  level  and  v"  refers  to  the  lower  level. 

Of  course  the  sum  of  iK)tential  energy  of  the  nuclei  and  the  electronic  energy  of  tlie 
electrons  forms  the  potential  well  in  which  the  nuclei  carry  out  their  vibrations.  It  turns 
out  that  this  potential  well  is  harmonic  to  first  order  (very  small  amplitude  vibrations 
alx)ut  the  equilibrium  displacement  ixitween  the  nuclei)  but  is  anharmonic  if  larger 
displacements  from  equilibrium  arc  considered.  This  anharmonic  effect  can  l)e  taken  into 
account  by  using  a  potential  of  the  form  (instead  of  V  =  ^); 

v  =  ^-gx3,  (2.14) 

where  x  =  r-r^,  r^  is  the  (xiuilibrium  separation  of  the  nuclei,  and  g  is  much  smaller  than 
k/2.  If  this  potential  is  substitutetl  into  the  Schodinger  equation,  the  following 
anharmonic-oscillator  wave  numirers  (energies)  are  obtained: 

^  =  G(v)  =  w^(v+l/2)  -  1/2)2  +  i.,^Y^(v+l/2)3  +  ...  (2.15) 
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Previously,  we  considered  Uie  rotational  eiiergit's  of  a  molecule  on  the  basis  of  it 
being  a  rigid  rotator.  It  is  quite  obvious  that  if  the  molecule  is  vibrating  it  is  not  rigid.  In 
a  rather  involved  wave  mechanical  calculation,  Pauling  and  Wilson*  show  that  if  one  takes 
into  account  the  molecular  vibrations  or  rotations  the  following  wave  numbers  for  a  given 
rotational  (J)  and  vibrational  (v)  level  are  given: 


l.\,(j)  =  H,.J(.i+])  -D,J2(J+1)^ 

(2.JG) 

where  Bv  =  Pe  “  Oe(v+l/‘2), 

(2.17) 

Dv  =  De  +  /4(v+l/2), 

(2.18) 

lie  =  (:j  (sex’ Fxpiatioii  2. 10), 

(2.10) 

1)  - 

(2.20) 

Oe  and  /ie  are  much  smaller  than  H*.  and  l)e  respectively. 


Finally,  we  consider  tlie  electronic  energy  states  as.sociated  with  the  electrons  about 
the  nuclei.  'I'his  energy  is  obtained  by  considering  the  electrons  interacting  with  each  other 
and  the  nuclei.  'I'liere  is  no  exact  solution  of  the  Schbdinger  e(|uation  for  these.  One  can 
('siimate  the  energies  of  each  level  for  a  given  molecule  using  the  Porn-Oppeidicimer 
aj)])roximation.3  However,  in  practice  the  dwtronic.  energies  are  observed  s{K;ctroscopically 
and  used  to  give  the  base  energy  of  each  electronic  level.  In  terms  of  wave  numbers  this 
energy  is  given  by: 

d’e  =  T„  +  AAV,  (2,21) 

where  T„  repre.sents  the  energy  of  the  electronic  state  (base  energy)  nc'glecting  energy 
splitting  of  levels  due  to  electron  spin,  A  is  the  spin  orbit  coupling  constant,  A  represents 
the  component  of  electronic  orbital  angular  momentum  along  the  internuclear  axis,  and  ^  is 
the  sum  of  electron  spins  for  the  molecule. 

If  the  electronic  orbital  angular  motnentum  is  giv('n  by  a  vector  L  with  magnitude 
L=|L|  then  the  allowed  quantum  nmniK'rs  representing  the  component  of  electronic 
angular  momentum  along  the  internuclear  axis  are  given  by: 
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A  =  0,1,2, . L  (2.22) 

One  should  note  that  A  =  0,  1,  2,  ...  corresponds  to  the  molecular  state  designation  E,  n, 
A,  ...  (e.g.,  NS(B'^n)  refers  to  nitrogen  sulfide  in  a  state  with  A  =  1).  'I'o  in  equation  (2.21) 
above  is  also  a  function  of  A. 

If  the  total  electron  spin  associated  with  the  molecule  is  S  and  the  magnitude  of,  this 
spin  is  given  by  S  =  |S),  then  S  will  |)reccs.s  about  the  fi(*ld  direction  (in  this  case  the 
interuuclcar  axis)  with  a  constant  comix)ncnt  E(h/2;r).  Quantum  theory  allows  the 
following  values  of  S  =  S,  S-I,  ...  -S.  (2,23) 

The  total  number  of  E  values  for  a  given  S  from  Equation  (2.23)  is  given  by: 

M  =  2S+1,  (2.24) 

where  M  is  called  the  multiplicity  of  the  state.  The  mnltiplet  structure  observed  in 
electronic  bands  is  due  to  the  electron  spin  and  the  slightly  different  energy  levels 
corresponding  to  the  different  eUx^tron  spins  given  in  Equation  (2.23).  'I'he  electronic 
energy  of  the  multiplet  is  given  by  Eijuation  (2.21).  Note  that  if  A  =  0  for  a  E  state,  there 
can  be  no  spin  related  splitting  or  multiplet  structures.  'I'his  is  bocau.se  when  A  =  0  there 
is  zero  orbital  angular  momentum  of  the  electrons  and  therefore  no  axial  magnetic  Reid. 
Since  this  field  causes  the  spin  splitting  there  can  be  no  splitting  even  if  there  is  a  non-zero 
multiplicity.  So,  a  state  luis  a  multiplicity  of  M=2,  but  because  A=0  there  is  no 
multiplet  splitting. 

The  comjwnent  of  the  total  angular  momentum  of  the  electrons  along  the 
internuclear  axis  is  given  by: 

n=|A  +  E|.  (2.25) 

This  plus  the  quantum  number  for  rotation  (call  it  N)  of  the  nuclei  gives  the  resultant  total 
angular  momentum  J.  If  the  interaction  of  the  electronic  motion  and  the  nuclear  rotation 
is  small,  but  the  electronic  motion  is  strongly  coupled  to  the  line  joining  the  nuclei,  0  is 
well  defined.  This  is  typically  referred  to  as  a  lluiid's  case  a.  There  arc  a 

() 


variety  of  eases  (see  llertzbcrg);  however,  for  NS  Hand's  cas(  a  ai)plies.  So,  combining  all 


the  energy  (wave  number)  terms  for  a  molecule  we  obtain  the  following  sum: 

=  Te  +  G(V)  +  F(J)  (2.26) 

The  wave  numl)Oi  associated  with  the  emission  of  a  photon  or  the  transition  from  one  state 
(n)  to  another  (m)  is  given  by; 

=  'I'e'  -  Te"  +  C(v')  -  C(v")  +  F(J')  -  l'(.J") 

From  Fapiations  (2.15-2.20),  this  beconurs: 

=  Te'-Te"  +  +  l v'' +  1 /2)  +  a»e'' ye''(v''+ 1 /2)H... 

+  Hv'J(.i'+i)  -  i)v'.)'‘^(.i'+i)M}v".)"(.)"+i)+i)v".r'^(.r'+i)^+....  (2.27) 

This  corresponds  to  the  wave  number  of  the  emitted  plioton  if  tlie  transition  is  allowed. 
Using  (luantum  mechanics  it  is  straightforward  to  show  that  molecular  transitions  from  one 
particular  state  to  another  are  generally  allowed  if; 


AJ  =  J'-J"  =  0,  ±1 
and  for  Hand's  case  a: 

AA  =  A'-A"  =  0,  ±1 

A^  =  =  0. 


(2.28) 


(2.20) 

(2.80) 


Note  that  Av  =  v'-v"  is  not  |■('st^icte(l  in  any  way. 

So,  if  we  apply  this  to  the  i)r()i)oscd  NS(H‘'^H)  -  NS(X2H)  transition,  we  first  note  that 
since  the  upixrr  and  lower  states  are  ^H's,  A=1  for  both  and  AA=0.  Since 
M=multiplicity=2S+l=2  =»  S=l/2  or  X'=±l/2.  So,  12=  |  A+X|  =  l/2,3/2.  However,  since 
A^=0  in  Hand's  cruse  a  the  possible  electronic  transitions  are  restricted  to  the  following; 
y  -  X2Hi/2  and  U2H3/2  -  X2H3/2, 

where  12  has  been  written  as  the  subscript  on  the  H  (orbital  angular  momentum 
corresponding  to  A=l). 

Keeping  in  mind  that  AJ=0,±1  in  Equation  (2.27),  one  can  represent  the  possible 
transitions  for  the  NS(B2ri)-NS(X2H)  graphically  (sa;  Fig.  1).  Hi,  Qi,  Pi  corresi)ond  to 
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the  branches  with  AJs  of  +1,  0,  -1  respectively  for  the  transitions.  R2,  Q2,  P2 

correspond  to  AJs  of  +1,  0,  -1  respectively  for  the  ^113/2^113/2  transitions. 


11.2.  POPULATION  OF  ENFllGY  LEVELS 


Up  to  this  point  we  have  consi<l('r(!d  the  electronic,  rotational  and  vibrational  energy 
levels  available  to  a  molecule.  We  have  also  considered  the  allowed  transitions  from  one 
state  to  another.  Since  we  have  determined  the  energy  levels  available  wo  now  need  to 
figure  how  many  molecules  populate  any  given  level.  Often  in  nature,  systems  exist  in 
states  of  tliermal  C(|uilibrium.  I’liis  is  typically  true  of  the  molecules  we  study.  For  this 
reason  we  first  consider  obtaining  the  populations  of  electronic-ro-vibrational  energy  levels 
for  molecules  in  thermal  equilibrium. 

Electronic  Distribution  Function 

First,  we  consider  the  electronic  distribution  function  for  a  doublet  state  in  thermal 
cxjuilibrium  (we  consider  doublets  because  of  the  doublet  NS  transitions).  For  a  doublet 
state  in  equilibrium  the  number  of  molecules  in  the  upper  electronic  state  with  a  given  il  is 
given  by  (we  are  interested  in  the  ui)per  state  because  we  want  to  produce  an  emission 
spectrum)  a  Uoltzmami  distribution; 


-6 


Ne 


(2.31) 

where  is  the  energy  in  wave  numbers  of  the  different  electronic  levels  associated  with 
different  values  of  il  (with  6min  =  0)-  Miu  is  the  minimum  f)  value  and  k  =  .695 
crn'*/“kelviii.  For  NS  we  have 

Ne-ei/2/kTj,  ^ 


N 


Ee  ^^/kTp^  e"^‘/2/‘^'‘'E+e 


(2.32) 


and  with  Gmiti  =  G1/3  =  0,  G3/2  =  90,  k'fp  =  200  ciir'  at  room  temperature,  where  Tp  = 
300ok  =  the  electronic  temperature. 
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Ne,/3=N(.61) 


N 


E3/2 


=  3/2 /kTp,  ^  N(.3893) 

So,  for  NS  in  the  B  state  the  0=1/2  level  has  5/3  the  i)opulation  of  the  0=3/2  level. 
Vibrational  Distribution  Function 

For  vibrational  levels  in  thermal  equilibrium  the  population  of  the  different  energy 
levels  is  given  very  accurately  by  a  Boltzmann  distribution.  The  number  of  molecules  in  a 
given  vibrational  level  is 

=  (2.33) 

where  Qv  =  ^  )hc/kTv^  ^  ^2.34) 

=  state  sum, 

Go(v)  =  Cl(v)-G(0),  (2.35) 

Tv  =  the  vibrational  temperature. 

It  should  be  noted  that  if  we  are  concerned  with  emission,  Go(v)  and  Nv  correspond  to  the 
upper  state.  The  G(v)  above  is  given  by  Equation  (2.15). 

Rotational  Distribution  Function 

The  number  of  molecules  found  in  a  given  rotational  energy  level  in  thermal 
equilibrium  is  given  by  a  Boltzmann  distribution  function  times  the  degeneracy  of  the  level 
which  is  usually  2J+1.  This  is  because,  for  a  given  energy  each  of  2J+1  degenerate  levels 
may  be  populated.  The  number  of  molecules  in  a  given  J  level  is  given  by; 


-F(J)hc/kTr 


where  F(J)  is  given  by  Equation  (2.16), 

Q,  =  1  +  +  ... , 

Tr  =  the  rotational  temperature. 

Now  F(J),  to  good  approximation,  can  often  be  set  cxpial  to: 
F(J)  =  Bv  J(J+1). 


(2.36) 

(2.37) 
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If  this  is  substituted  into  (2.36)  and  (2.37), 

Q,  =  1  +  3(,-2Bvhc/kTr  ^  ^^-(iUvlic/kTr  ^  ^2.38) 

k'F 

“  ITcl^  ^  Bv,  which  is  usually  the  case.  (2.39) 

Tr  is  usually  greater  than  300oK  and  IK-  is  about  .6  for  NS.  This  makes  (2.39)  a  good 


approximation.  This  implies  that 

\r  ~  NhcBy  (2.1  +  1  '  )hc/k'rr 


(2.40) 


At  this  i)oint  we  are  in  a  position  to  fairly  accurately  define  Nn  of  Equation  (2.8)  by  using 
(2.31),  (2.33)  and  (2.36): 


n'  of  course  corresponds  to  the  upiKtr  state  with  fi',  J',  v'. 


(2.41) 


II.3.  ELECTRONIC-RO-VIBRATIONAL  INTENSITY  (COUNT)  FUNC'J'ION 

We  almost  have  everything  wo  need  to  define  E(iuation  (2.8)  for  the  NS  molecule. 
Nr,  of  (2.8)  is  given  by  Equation  (2.41),  of  (2.8)  is  given  by  (2.27).  d,,  of  (2.8)  is  C(|ual 


to  the  degeneracy  of  the  levels  which  is  2.1+ 1.  In  the  literature 


5:|R"i'”kr'  =  iRv'v'T'  Mv'v"  Sj,  j.„  (2.42) 

where  the  R^'v  is  the  part  of  the  transition  probability  depending  on  the  electronic  and 
vibrational  eigenfunctions  and  Sjijk  is  the  part  of  the  line  strength  that  depends  on  .1  (the 
total  angular  momentum)  and  is  called  the  Honl-bondon  factor.^  In  the  following  we  will 

show  how  Rv'v"  ~  qv'v"  and  define  all  terms.  'I'he  electronic  and  vibrational  part  of 


the  eigenfunction  can  be  written  as: 


X  =  XeXv 


(2.43) 


The  probability  of  a  transition  betwerni  two  different  (as  shown  in  Eriuation  (2.2'-4)) 
electronic-vibrational  states  is  given  by: 
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|R|2=  I  f  ;t'Mxdr|2  (2.44) 

The  only  parts  of  the  electric  moment  M  that  can  contribute  to  R  are  those  that  depend  on 
the  coordinates  of  the  electrons.  Rewriting  (2.44)  and  substituting  (2.43)  with  this  in  mind 
gives: 

R  =  JMc  JCe'JTv'XeK*'’' 

where  Xy  correspond  to  the  vibrations  of  the  nuclei  along  r  and  therefore  only 
contribute  along  r,  and  Xq  de{)end  on  the  electronic  coordinates  only  and  contribute 
accordingly. 

The  first  integral 

ft.  =  /  A:e''MeV''’'c  ' 

is  called  the  electronic  transition  moment.  This  transition  moment  deiKUids  only  slightly 
on  the  internuclear  separation  of  the  nuclei.  For  this  reason  Rg  is  often  ap])roximated  by 
its  average  value  over  a  transition: 

Re  =:  Re  (2.47)* 

The  second  integral  of  (2.45)  is  called  the  overlap  integral.  It  basically  states  that  a 
transition  is  more  probable  if  the  wave  functions  simultaneously  have  large  values  for  the 
two  vibrational  states  at  a  given  internuclear  .separation.  This  embodies  the 
Franck-Condon  principle.  As  a  result  the  overlap  integral  is  often  written  in  this  way: 

where  Qv'v"  is  the  Franck-Condon  factor. 

Note,  .sometimes  one  must  include  the  variation  He(i')  with  f  to  accurately  model  the 
spectrum.  This  happens  to  be  the  case  for  NS.  As  a  function  of  the  r-centroids  (f^iyii)  of 
the  molecule  Re  is  given  by:^ 

Rc  =  lRe(f)|  =  Const.(-l+1.767rv'v"  -  1.037f2v'v"  +  .202  r3v'v")  (2.49) 

where  rv',v"  is  the  average  separation  of  the  nuclei  for  the  v'  to  v"  transition. 
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Now,  all  the  tertiis  are  defined  necessary  to  give  the  intefisity  of  a  spectral  line  for  an 
electronic-ro-vibrational  transition  from  one  state  to  anotlier  in  equilibrium.  Substituting 
(2.41,  42)  into  2.8  wo  liave: 

term  N  term 

0  V 

jJ'v"  _G4r'N 

‘j"v"n  -'W~- 

N  I  term 

13  M'  +  l)hc/kTr' 

.  //Jl{e^qv'v"S  j,  jM  (2.50) 

l-'or  NS,  we  first  used  lC(|uation  (2.50)  to  obtain  the  e(|uilibrium  synthetic  spectrum  at  high 
pressures.  The  synthetic  si)ectrum  closely  matched  the  experimental  si)ectrum. 

The  low  pressure  cases  (<2  'lorr)  had  nou-e(iuilil)rium  rotational  and  vibrational 
distributions.  We  substituted  our  guesses  for  the  distributions  (Ny  -f-  Nj  terms  in  (2.50) 
into  the  code  until  the  synthetic  spectrum  matched  the  experimental  spectrum.  In  this 
way  we  determined  the  non-c(|uilil)rium  vibrational  and  rotational  population  distributions 
of  the  experimental  spectrum.  Such  a  comparison  is  shown  in  Figure  2,  in  which  the 
relative  populations  of  the  vibrational  levels  v'  =  0-G  (for  n21l|/>  -  X'^lli/o  and  - 

X2n3/2  transitions)  are  given.  The  Nj  term  is  given  in  FUNCTION  XI  of  the  code,  run  at 
a  rotational  temperature  of  1400"K. 

III.  HOW  TO  RUN  THE  CODE 

To  run  the  NS(B2n)  -  NSfX^H)  spectrum  code  one  needs  to  first  get  in  the  directory 
where  it  resides.  It  is  stored  there  along  with  files  filled  with  Franck-Condon  factors 
(FRANK. CON),  Dunham  coefficients  (DUNHAM.(X)F),  r^entroids 
(RCENTROID.DA'I’),  the  population  weights  (POIM/I'N.WT),  file  of  data  file  names 
(FILE.OMA),  and  the  experimental  data  files  (named  in  FILE.OMA).  These  files  must  be 
present  when  running  the  program. 


e'^D/k  r,,  [o'^Jo(v' )hc/kT;' 

^  • - - 

Ee  ^f/k'l'.J 

a  I  j 


First  one  must  login  to  the  VAX  as  follows: 

C  FJ 

USER:  LAWCONNELL 

PASSWORD;  ***  ,  see  Captain  Lawconnell. 

Then  default  the  directory  to  ASPECTRUM  by  typing: 

Sd.ASPECTRUM 

Assuming  all  the  Dunham  coefficients,  Franck-Condon  factors,  etc.,  mentioned  above  are 
available,  all  one  need  do  is  type: 

RUN  SPNS 

This  activates  the  executable  code  found  in  file  SPNS.EXE.  The  program  is  interactive 
and  the  input  required  is  self-explanatory.  The  user  will  be  asked  to  input  the  device  on 
which  to  plot  the  spectrum,  the  spectrum  title,  the  rotational  and  vibrational  temperature, 
whether  it  is  desired  to  plot  experimental  data  versus  the  theory,  etc.  If  one  wishes  to 
input  non-equilibrium  population  weights  for  the  upper  state  vibrational  levels,  this  is 
done  by  modifying  the  file  POPLTN.WT  and  specifying  that  a  non-equilibrium 
distribution  is  required  when  asked  while  doing  the  input. 

There  is  a  file  called  INPUT.DAT  that  saves  all  the  input  parameters.  If  one  is 
running  a  variety  of  cases  that  are  similar  this  file  will  be  accessed  to  give  the  input.  The 
code  prompts  for  any  new  changes.  In  addition  the  code  produces  three  output  files;  1) 
CHECK-OUT  which  contains  the  vibrational  transition  lines  -  and  other  information  that 
indicates  whether  or  not  all  is  well  with  the  run,  2)  OUTPUT.SPEC  which  contains  a  list 
of  intensities  versus  wavelength  and  3)  the  plot  file  (as  sixicified  by  the  user). 

Most  features  of  the  code  are  documented  internally.  The  file  structure  for  the 
various  data  files  mentioned  earlier  can  easily  be  determined  by  examining  the  respective 
read  statements  in  the  code.  This  also  applies  to  the  output  files. 
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If  one  needs  to  transfer  data  from  the  OMA  to  the  VAX  one  way  to  do  it  is  to  take  a 
floppy  of  the  data  (DOS  format)  and  use  a  PC  along  with  KHHMIT  to  transfer  the  data. 
The  procedure  for  this  is  outlined  Ix'low: 

Type: 

ZS 

KERMIT 

SERVER 

Hit  the  <ALT>  key 
Typo; 

K 

S  (to  send) 

A;*.* 

*  * 

F  (to  finish) 

If  the  procedure  failed  start  over  by  typing: 

<AI/r> 

Exrr 

VAX 

EXIT 

and  then  start  with  the  KERMIT  command  and  procctnl  as  before. 

Finally,  if  one  needs  to  modify  the  code  to  do  some  problem  of  interest  it  is  necessary 
to  know  how  the  previously  presented  c(iuations  relate  to  the  code.  If  we  start  at  the  top, 
the  first  routine  encountered  is  the  program  driver.  In  it  are  documented  all  the  references 
used  in  writing  the  code  and  for  producing  the  input  files.  The  program  driver  is 
responsible  for  the  program  flow  from  plot  driver  initialization,  through  the  six)ctrum 
generation,  to  the  output  of  the  siK’ctrum. 


If) 


The  next  routine  is  called  SUBROUTINE  INPUT.  It  reads  in  the  experimentally 
derived  Dunham  coefficients.  Dunham's  equation  combines  Equations  (2.15-21)  into  one 
equation; 

T  =  Y/j(v+ 1  /2)j  jj( J+ 1  )j  (2.51 ) 

where  the  Y^j  are  the  Dunham  coefficients.  From  Equation  (2.15-21)  it  is  easily 
determined  that  the  Dunham  coefficients  correspond  to: 

I  Yoo  =  Te  Yoi  =  Be  j  Y02  =  ~De 

^  Yjo  =  We  Yii  =  “Qfe  I 

I  ' 

^  Y20  =  -WeXe  •••  I 

These  arc  the  values  found  in  the  file  DUNHAM. COF.  We  typically  only  use  the  values  in 
the  dashed  box  to  model  the  NS  si)cctra  to  the  resolution  of  our  spectrometer.  It  should  be 
noted  that  there  is  a  different  group  of  Dunham  coefficients  for  each  electronic  level  -  four 
groups  for  NS.  This  is  referenced^  and  explained  in  greater  detail  within  the  code. 

The  Franck-Condon  factors  ((iv'v")  of  Equations  (2.48)  and  (2.50)  are  given  in  file 
FRANK.CON.5  These  aie  read  in  SUBROUTINE  INPUT.  In  addition  the  f-centroids^ 
(fv'v")  of  Equation  (2.49)  are  also  road. 

If  a  non-equilibrium  vibrational  distribution  is  needed  the  Ny  term  of  Equation 
(2.50)  is  replaced  by  the  population  weights  found  in  file  POPLN.WT.  These  are  also  read 
in  INPUT.  In  addition  the  rotational  temperature  (Tr)  and  vibrational  temperature  (if 
equilibrium)  are  read  in  for  the  Nj  and  Ny  terms  of  (2.50)  respectively. 

Go(v),  Qy,  By,  Tg,  (Equations  2.15,  35;  2.34;  2.17;  2.32;  2.32  respectively)  are 
obtained  in  SUBROUTINE  INTENSE.  These  values  are  used  there  to  calculate  the  Ng, 
Ny  and  Nj  terms  of  Equation  (2.50).  All  but  the  rotational  part  of  Equation  (2.27)  is 
calculated  in  this  routine  and  stored  in  matrix  TRW. 
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The  rotational  part  of  (2.27),  taking  into  account  the  P,  Q,  and  R  branching,  is  done 
in  SUBROUTINE  DOUBLET.  This  allows  for  the  determination  of  the  i/  in  equation 
(2.50).  DOUBLET  then  calls  SUBROUTINE  AINT  which  gives  the  Ildnl-London  factors 
(Sj'j")  of  equation  (2.50). 

Since  the  final  intensity  is  normalized,  the  actual  value  of  N  and  the  constants  of 
equation  (2.50)  are  not  important.  All  the  terms  mentioned  in  (2.50)  arc  brought  together 
in  FUNCTION  XL  This  is  where  the  intensity  of  a  given  spectral  line  is  calculated.  These 
intensities  are  redistributed  according  to  a  Lorentzian  slit  width  broadening  redistribution 
function  in  SUBROU'IMNE  LORN'l'ZN: 


where  7  =  the  spectral  resolution  in  Angstroms  (which  is  on  the  order  of  two  .\ngstroms 
for  our  OMA  so  that  Doppler  broadening  was  not  considered), 

A  =  the  wavelength, 

I  X  10®/i/j'j"  =  A  j,  =  the  wavelength  of  the  given  transition  one  is  redistributing, 
lo  =  original  intensity  of  the  line. 

The  experimental  data  is  read  in,  spliced  and  weighted  in  SUBROUTINES 
LOADFIL,  SPLICE  and  WEIGHT,  respectively.  The  data  and  the  synthetic  spectra  are 
normalized  and  plotted  in  SUBROU'PINE  OUTPUT  (as  specified  by  the  user). 


IV.  NS(B  -  X)  SPEC'l’RUM  CODE  LISTING,  FILES  AND  SAMPLE  SPECTRUM 


What  has  teen  written  above  is  embodied  and  takes  a  functional  form  in  the 


following  program  listing.  In  addition,  the  files  used  to  run  the  coile  arc  also  included 
(they  follow  the  code).  Finally,  a  sample  spectrum  (Fig.  2)  comparing  the  code  output  to 
experiment  is  included.  It  should  be  noted  that  the  theory  and  exi)eriment  correspond 
fairly  well.  This  was  true  of  all  the  s|X!ctra  generated  with  the  code  and  compared  to 
experiment  (for  a  variety  of  rotational  tem|X}ratures  and  pressures). 
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PPOGR/i»»  SfECTSf 


WRITTEN  EY: 

CAPT  eOB  LAWCONNELL 

FRANK  J.  SEILER  RESEARCH  LABCRATORY 
UNITEC  STATES  AIR  FORCE  ACAOEPY 
COLORADO  SPRINGS^  COLCRACC  SC84C 
PHONE:  719-472-35C2/  AV  259-35C2 

THIS  PRCGRAP  COPFUTES  ANC  PLOTS  THE  NS(6)  TC  NS()(>  SPECTRUM 


REFERENCES:  1)  G.  HERZBERG/  SPECTRA  OF  CIATOMC  MOLECULES^  VNR  C  19SC^ 
F. 201/268-271 . 

-  RELATION  GIVING  INTENSITY  OF  LINE  FCR  TRANSITION 

2)  P.  BCRN  AND  R.  OPPENHEIHER/  ANN  PHYSIK  84/457(1927). 

-  aCRK-OPPENHEIMER  APPROXIMATION 

3)  OBASE/ET  AL/  J  CNEM  PHYS  89/257(1988). 

-  FRANCK-CONOCN  FACTORS/  R  CENTROIDS/  ETC.  FOR  NS 

4)  J.  JEFFRIES/  ET  AL/  J  PHYS  CHEM/  IN  PRESS/  1  DEC  87 

-  EXPERIMENTALLY  MEASURED  TRANSITION  PROBABILITIES 
(IE.  R£**2). 

5)  R.  EHGLEPAN  Jfi./ET  AL/  LA-4364/UC- 34/1 ID-4500. 

-  AIM  SUBROUTINE  FOR  HONL-LCNDON  FACTORS  IN 
DCUELET  TO  COUBLET  TRANSITIONS 

6)  I.  KCVACS/  ROTATIONAL  STRUCTURE  IN  THE  SPECTRA  OF 
OIATCPIC  MOLECULES#  APER.  ELSEVIER  PUE/  OLT  OF  PRINT/ 
1969. 

-  BASIS  CF  AIM  SUBROUTINE 

7)  K.  P.  HUEER  AND  G.  HERZBER6/  MOLECULAR  SPECTRA  ANC 
MOLECULAR  STRUCTURE/  VOL  IV/  VNR  C  1979/  P486-487. 

-  DUNHAM  COEFFICIENTS  FOR  NS 
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PARAPETGR(NRCT=125) 

THE  COMMONS  ARE  DIMENSIONED  AS  FOLLOWS  (THEY  MUST  EE 
DIMENSICNEC  GREATER  THAN  OR  EQUAL  TO  THE  INDICATED 
VARIABLES): 

CCMMON/FACT/F(NF/20/20)/T(ND/NOl/2)/G(ND/2C)/ 

1  e(ND/20)/DS(NC/2C) 

CCMMON/P  /TRW  (6/20/70 /FI  (6/NROT  ) /P2  (6/NPOT)  / 

1  FI  INI  (6/NfiCT) /P2lNT(6/NRCT) 

COMMON/Q  /Q1 (6/NRCT) /Q2(6/NfiCT)/ 

1  Cl INT(6/NR0T)/a2INT(6/NR0T) 

CCMMON/R  /R1 (6/NROT) /R2<6/NRCT)/R1INT(6/NR0T)/ 

1  R2 INT (6/NROT) / FIN T(6/1COCC)/FIN(1COOC)/FINX(1COOC) 
CCMMON/TEMP/TV/TR 

WHERE/  NF«NT*THE  NUMBER  OF  GROUPS  CF  FRANK-CONOON  FACTORS  -- 
1  6RCUP  PER  ELECTRONIC  TR A NSI T ION/ 1 . E .  B3/2  TC  X2/2 
AND  el/2  TO  X1/2  FOR  NS. 

NF1»  V  PRIME  MAX  INDEX  (UPPER  STATE) 

NF2>  V  DOUBLE  PRIME  MAX  INDEX  (LOWER  STATE) 

ND»  NUMBER  CF  DUNHAM  COEFFICIENT  GROUPS  -- 
ONE  FOR  EACH  OMEGA  6UANTUM  NUMBER 
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C  N01»  THE  L  VALUE  F0«  THE  OUHHAH  COEFFIEKT  HATRI*: 

C  Y(L^J) 

C  LWAV=  MINIHL'F  WAVELENGTH  CF  SPECTRM  <AN6STRCHS> 

C  LWAV1=  HAXIFLM  WAVELENGTH  OF  SPECTRUH  (ANGSTROMS) 

C  NR0T»  NIPBER  OF  POSSIBLE  ROTATIONAL  LEVELS 

C  NOTE:  ALL  INTEGERS  USEC  TC  DEFINE  ARRAYS  REFERIN6  TO 

C  VIBRATIONAL  LEVELS  ARE  ONE  LARGER  THAN  THE 

C  CORRESPCNOIKG  VIE.  LEVEL. 

C 

CHARACTER*A5  RAHCEV 
C 

CALL  ASSIGN<3/*CH£CK.0UT*) 


WRITE(E /*) 

*UHA 

T  TYPE  Of 

DEVICE 

DO  YOU  WANT 

TO  PLOT 

URITE(€#*) 

•SPE 

CIFY:  1 

FCR 

HP  755C* 

URITE(E^*) 

• 

2 

FOR 

TEKTRONIX 

4107* 

WRITE(6^*) 

1 

2 

FCR 

VT240* 

URITE(E^*) 

• 

4 

FOR 

VT  IOC  RETRO* 

URITECE/*) 

« 

5 

FCR 

TEKTRONIX 

401C* 

WRITECE/*) 

t 

E 

FCR 

ANY  ASCII 

PRINTER* 

REA0(6#«)  I 

PLOT 

IFdPLCT.EC 

.1) 

NAMOEV»*HP 

7550* 

IFdPLCT-EC 

.2) 

NAHCEVsMFKTRONIX 

41C7* 

IFdPLCT.E: 

.3) 

NAMCEV** VT 

240* 

IFdPLCT.EC 

.4) 

NAMCeV**VT 

ICO  RE 

TRC* 

IFdPLCT.EC 

.5) 

NAMCEV**TEKTRONIX 

4C10* 

IFdPLCT.EC 

.6) 

NAHCcV»*ANY  ASCII 

PRINTER* 

CALL  PLTOEV 

(NAP 

CEV) 

CALL  INPUT 
CALL  INTENSE 
CALL  OL'TPUT 
CALL  DCNEPL 
CALL  CLOSE(3) 

STOP 

END 

SUBROUTINE  INPUT 
PARAHETER(NROTa125) 

C0PHCN/FACT/F<3/3Q^30)#T(9^3C^2)#G<9/3C)/ 

1  3<9^3C)#0G<9/3C) 

COPHCN/P  /TRW <6# 30^ 30)^ PI (6# NROT) / P2 (6# NRCT) / 

1  P1INT<6^NRCT)^F2INT(6^NR0T> 

COPHCN/G  /Q1(E/NR0T>#Q2(E#NRCT)^ 

1  61INT(6^NROT)^62INTC6#NROT) 

COPHON/R  /R1(e^NR0T)#R2(E#NRCT)^RlINT(6#KR0T), 

1  R2INT(6/NROT)^FINT(6^tCOOC>^FIN(1CCOO)/FIHX(1CCOO) 
COPMCN/TEMP/TV/TR 

COPPCN/INPLT/NTITLE^IFDIS^ieRANCH#Wr:N#WPAX/RRE<^ 

1  XTV^XTR^IFDCNE^IFCEF^IPIN^IPaX^IENVH^ENVHT^V# 

2  XP^TK/GANPA2 

C0PN0N/1NDEX/NF^NF1^NF2^NT^ND^ND1 
C0PP0N/U6HT/NGHT  (3C/'30)^NVIE^IFEXCIT 
COPNON/CNTR/CNTR  (3/ 30^30  #1F CENT  R 

COPHCN/EXPOAT/IFDATA^WAVE(2C^1C30)^SPEC(20/103C)/UAV(2CCOO)^ 

1  SP(20CCO}^WW(2COCC}^SPWC2CCCO)^WDAT(2C)^SEN(2C)#ICCUNT/IIPAX# 

2  INNAX/BKGRND/SLOPE 

C0RHCN/PERT/IFPCRT/TRP<IFPLB^VAR^VAR1^VAR2 
CHARACTER*AC  NTITLE 
CHARACTER*3C  IFILE 
C 

c 

C  READ  ZN  THE  APPROPRIATE  FRANK  CONDON  FACTORS 
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c 

CALL  ASSIGN (2^* FRANK. COK*  ) 

READ(2/*)  NF^NF1^NF2^NT 
WRITE(2^*)  NF^NF1i-NF2/-NT 
DO  10  Nsl/NF 

C  N=1  inPLIES  81/2  —  X1/2  ELECTRCNIC  TRANSITION 
C  N=2  IMPLIES  83/2  ~  X3/2  ELECTRCNIC  TRANSITION 
DO  20  I  =  UNF1 

READ(2^*)  (P(N,l,If),IP=1,nf2) 

WR1TE(3#*)  (F(N,I^1P)/IP=1^NF2) 

20  CONTINUE 
10  CONTINUE 
CLCSE(2) 

C 

C  READ  IN  THE  DUNHAM  COEFFICIENTS 
C 

C  NOTE:  THE  ZERO  POINT  ENERGIES  ARE  INCLUDED  IN  THE  FIRST 
C  DUNHAM  COEFFICIENT  IN  THIS  FORMULATION 
C 

CALL  ASSIGN(2^*CUNhAM.CCF*) 

READ(2^*)  ND/NDI 
WRITE{3/*)  N0#>NC1 
DO  30  ►sl^NC 

C  M=1  IMPLIES  X2FI1/2  STATE 

C  M=2  implies  B2F11/2  STATE 

C  M=3  IMPLIES  X2FI3/2  STATE 

C  M»<  IMPLIES  B2PI3/2  STATE 

DO  40  L=1/N01 

REA0(2/*)  Y<M/L^1)/Y<H/L#2) 

40  CONTINUE 

DO  39  L=1/N01 

WRITE<3^a>  Y<M/L^1)^Y<M/L/2) 

39  CONTINUE 
30  CONTINUE 

CALL  CLCS£(2) 

WRITE<E<*)  'IF  YCU  HAVE  NOT  DONE  THIS  SPECTRUM  BEFORE  CR  WOULD*/ 

1  '  LIKE  TO  START  FROM  SCRATCH  TYPE:  C  /  OTHERWISE  TYPE:  1  .* 
READ(6/*)  IFDQNE 

IFdFDCNE.EG.O)  THEN 

WRITECE/*) *WNAT  DC  YOU  WANT  TC  ENTITLE  YOIR  SPECTRUM?* 

READ(d/33)  NTITLE 
33  FCRMATCA3C) 

HRITE<6/*)  'IF  YOU  WANT  THE  CCDE  TC  DETERMINE  THE  MINIMUM  AND*# 

2  *  maximum  wavelengths  between  WHICH  THE  SPECTRUM  IS  TC  BE*# 

3  *  PLCTTEC  TYPE:  C  '# 

4  '  NOTE:  THIS  RESCLVES  THE  SPECTRUM  ON  20CC  POINTS  IN  WAVELENGTH* 
READCd#*)  IFDEF 

IF(IFOEF.NE.C)  THEN 

HRITE(6#*)  'WHAT  IS  THE  MINIMUM  W A V EL E NG T H ( WM IN)  IN  ANGSTROM?* 
REAC(6#«)  WMIN 

WRITE(6#*)  'WHAT  IS  THE  MAXIMUM  W AVELENGT HCWM AX)  IN  ANGSTROM?* 
REAC(S#«)  WMAX 

WRIT£(6#*)  'WHAT  RESCLUT ICNIRE  S)  DO  YOU  DESIRE  IN  THE  SPECTF??* 
WRITE<6#*)  'NOTE:<WMAX-WMIN)/RES  .LE.  10CO0  AND  MUST  BE  AN*# 

1  'INTEGER* 

REAC(6#«)  RRES 
ENDIF 

WRITE<6#*)  'WHAT  IS  THE  VIBRATIONAL  TEMPERATURE  (DEG  KELVIN)?* 
READCd#*)  XTV 

'  HRITE(6#*)  'WHAT  IS  THE  ROTATIONAL  TEMPERATURE  (DEG  KELVIN)?* 
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READ(6/*)  XTR 

WRITE(6/*)  *0C  YOU  WANT  A  LORENTZIAN  (O  CR  GAUSSIAN  (1)*/ 

1  •DISTRIBUTION  {GAUSSIAN  NOT  YET  IR P LENE NTE 0 ) ? * 

READ(£^*)  IFDIS 
IF(IFCIS.EQ.O)  THEN 

tiRITE(6#«)  *h»<AT  IS  6ANNA«*2  FOR  LORENTZIAN?  IF  C  IS  RETURNED*^ 
1  •GAHP A**2»9.56* 

READ(6/«)  GAFPA2 
IF<GAPPA2.EQ.C>  THEN 
GAMPAZsS.SS 
ENDIF 
ENOIF 

IF(  IFCIS.EC.1 )  THEN 

WRITE(6^*)  •WHAT  IS  THE  GAS  KINETIC  TEMPERATURE  (KELVIN)?* 
REAO<6^*)  TK 

WR1TE<6/*)  'what  is  the  EPITTING  POLECULES  HCUECOLAR  KEIGHT*# 

1  • (GRAPS/MCLE)?  • 

READ(6^*)  XN 
ENOIF 

WRITE(G^*>  'WHAT  IS  THE  EXPECTED  ENVELOPE  WIDTH  OF  TME*^ 

1  'THE  OISTRIBUTION  (IN  ANGSTROPS)--  IF  ZEfiC  IS  RETURNED  THE*^ 

2  'CODE  WILL  ESTIMATE  A  VALLE.* 

REAOCE^*)  ENVHl 

WRITE(G^*)  'WHAT  TYPE  CF  RCTATICNAL  BRANCHING  IS  REOUIREO?* 
WRITc(6^*)  *  THE  CHOICES  ARE:  * 

WRITE(6^*)  *  1)  Pl#Q1^fi1/F2/a2/R2  FOR  2PI--2PI  TRANSITIONS* 
WRITE(6/*)  '  CCRRESPONDING  TO  HLNDS  CASE  A.' 

READ  (G/*)  IBRANCH 

WRIT£(6/*)  *IF  YOUR  UPPER  ELECTRONIC  STATE  HAS  A  NCN-ECUIL I BRIL'M 
WRITE(6/*)  'VIERATIONAL  OITRIBUTICN  TYPE:  1  * 

WRITEC6/*)  '  OTHERWISE  TYPE:  C  * 

WRITE(6#*)  'HAVE  THE  RELATIVE  WEIGHTS  CF  EACH  VIBRATIONAL  LEVEL* 
URITE(6^*)  *IN  FILE  POPLTN.WT* 

REAO(E^*)  IFEXCIT 
IF(IFEXCIT.EQ-I)  THEN 
CALL  ASSI6N(2/*PCPLTN.HT*> 

REAO(2/*}  NVIE 

C  NVie  IS  THE  NUPBER  CF  VIBRATIONAL  LEVELS  IN  THE  UPPER  STATE  YOU  HAVE 
C  WEIGHTS  FOP.  ALSO  THERE  IS  ONE  SET  CF  WEIGHTS  FOR  EACH  TRANSITION. 

CO  SC  N-1#NF 

PEAO(2/>*)  (HGHT(N^I)^I*1#NVIE) 

WRITE(3/*}  (WGHT(NrI)#I«1^NVI6) 

SO  CONTINUE 

CALL  CLOSE(2) 

ENOIF 

WRITE(6/*)  *IF  YCU  WANT  TO  WEIGHT  THE  SPECTRAL  INTENSITIES* 
WRITE(6/*)  'WITH  THE  R-CENTRCIOS  TYPE;  1  * 

WRITE(6/*)  *  OTHERWISE  TYPE:  O' 

WRITE(6-f*)  *THE  R-CENTROICS  ARE  ASSUMED  TO  BE  IN  RCENTROIO.DAT* 
REAO(d^*)  IFCENTR 
IF(IFCENTR.EQ.I)  THEN 
CALL  ASSI6N(2/ 'RCENTRCID.OAT') 

DO  70  N«1/NF 

C  N*1  IMPLIES  BI/2  --  X1/2  ELECTRONIC  TRANSITION 
C  M«2  IMPLIES  B3/2  —  X3/2  ELECTRONIC  TRANSITION 
C  THERE  SHOULD  BE  AS  PANT  R-CENTRCIOS  AS  FRANK-CCNOCN  FACTORS 
00  80  I«1^NF1 

READ(2#«)  (CNTR (N/ I / IP) P« 1 /NF2 ) 

WRITECI/*)  (CNTfl(N/I#IP)#IP«1^NF2) 

80  CONTINUE 
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^0  COKTlNtE 

CALL  CL0SEC2) 

ENOIF 

WRITECE^*)  ‘INCLUDE  NON-0QLT 2MANN  C C NT R 1  BUT  I  ON  TC  SOTATIOKAL* 
WRITECE/*)  ‘LEVELS  AS  SPECIFIEC  BELOW  ?  YES  (1)/  NO  (0).‘ 
READ<6/*)  IFPERT 
IF (I FPERT.EQ.1)  THEN 

WR1T£{6^*)  ‘WFAT  IS  THE  FSEUCC  RCTATICNAL  TEHFERATURE  <K)?‘ 
READ(6^*)  TRF 

WRITE(6/’*)  ‘WHAT  IS  THE  RELAT  HEIGHT  CF  2N0  EKP  FOR  PCT  OIS?‘ 
FEAD(6/*)  VAR 

URITE(6^*)  'AT  WHICH  RCT.  LIKE  OCES  CCNT.  TAKE  EFFECT‘ 

PEAD(6^*)  VAfil 
VAR2=VAR1-1 
ENCl  F 

WRITE(6/-*)  ‘PUBLICATION  QUALITY  (GRAPHS  BLACK)?  YES  (T)/  NO  (0).* 
READ<6/*)  IFPUE 

WRITECdr*)  ‘DO  YOU  WANT  TO  PLOT  THE  THEORETICAL  ANC  E XP E R IP E NTA L * 
WftITEIdi.*)  ‘SPECTRA  ON  THE  SAME  GRAPH?  YES  (1)  /  NC  (C)‘ 

REA0(6/«)  IFOATA 
IF(IFOATA.EG.I)  THEN 

WRITE<6/*)  ‘SPECIFY  THE  THE  AVERAGE  NUPEER  OF  EACKGROUKC  CCL'NTS‘ 
WRITE(6,*)  ‘PER  PIXEL  TC  SLETRACT  FROP  YOUR  DATA.  CC  THIS  EY‘ 
WRITECEr*)  ‘FIRST  SPECIFYING  THE  INITIAL  NLP9EP  CF  COUNTS  TC‘ 
WRITE(C/*)  ‘SUETRACT  AT  YOUR  MIMKUP  WAVELENGTH* 

READ(6/*)  EKGRKC 

WRITECE/*)  'NCU  SPECIFY  THE  SLCPE  CF  THE  LINE  CF  THE* 

WRITECe^*)  ‘AVERAGE  INCREASE  IN  COUN TS / AN ST R OM ‘ 

REA0<6^*)  SLCPE 

WRITECE/*)  'YCLR  SPECTRAL  DATA  IS  ASSUMED  TC  BE  IN  THE  ‘ 
WRIT£<C/*)  ‘FILES  NAMED  IN  :  FILE. CPA.  NOTE  THAT  IT  IS* 
WRITECE/*)  ‘ALSC  ASSUMED  THAT  THE  FIRST  FILE  LISTED  IN  FILE.OMA* 
WRITECE/*)  ‘CORRESFCNOS  to  the  lowest  WAVELENGTH  ANC  SC  ON  ‘ 
WRITECE/*)  ‘IT  IS  also  ASSLPEO  THERE  ARE  NC  GAPS  IN  WAVELENGTH* 
WRITECE^*)  ‘BETWEEN  CONSECUTIVE  OMA  FILES.* 

ENOI  F 

IFdFDATA.KE.DTHEN 

WRITECi/-*)  ‘CALCULATING  THE  SYNTHETIC  SPECTRA...* 

WRITE (6^*)  ‘CRUNCH. . . C RUNC H . . . C RUN C H . . . ‘ 

ENDIF 

ELSE 

URITE<6r*)  ‘WHAT  IS  THE  NAPE  CF  YOUR  INPUT  FILE  (TYPICALLY 
1  * INPLT.DAT)?* 

READ(E^35)  IFILE 
CALL  ASSIGN(2/IFILE) 

READ(2^*)  WMIN/WMAX/RRES/XTV^XTR 
READ(2/35)  NTITLE 

35  FCRPAT(A3C) 

36  FCRPAT  (1  X/ A  30 

READ(2/*)  IFDI£rIERANCH#ENVH1/TK^XP^GAPPA2 
REA0(2/*)  IFEXCIT/ IFCENTR 
CALL  CL0SE(2) 

HRITE(6/*)  ‘THE  INPUT. CAT  FILE  HAS  THE  FOLLOWING  VARIABLE*/ 

1  ‘  VALLESi* 

WRITE<6/*)  ‘THE  PLOT  TITLE  IS  1:* 

WRITE(6/36)  NTITLE 

WRITE<6/*)  ‘LCRENT2IAN(0)  OR  GAUSSIANd)  CIS.  2:*/IF0IS 
WPITE(6/*)  'THE  ENVELOPE  HALF  WIDTH  CF  OIS.  !:*/ENVH1 
WRITE(6/*)  ‘THE  BRANCHING  IS  GIVEN  BY* 

WRITE(6/*)  *1)  2P1--2PI/HLNOS  CASE  A  A:*/IBRANCH 
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WRITE(6^*)  *THE  MINIMUM  URVELEK6TH ( A)  IS  5:*/UHIK 

WRXTE<6/*>  *TN£  IRAXINUR  WAVELENGTH  (A)  IS  6;*^WNAX 

IFOEFal 

WRXTE(6^*)  *THE  NIN  RESCLtiTXON  XN  ANGSTONS  7:*^RRES 
HRXTEEG#*)  *THE  ViBRATXONAt  TERP  (KELVIN)  IS  fir'/XTV 
WRXTE(6/*)  *THE  RCTATXONAL  TEPP  (KELVIN)  XS  9:*/XTR 
IF(XFDXS.EQ.O)  THEN 

HRXTE(6^*)  IS  1C;*/6AF!PA2 

ENOIF 

IFdFCXS.EQ.D  THEN 

WflXTE(6/*)  ‘THE  GAS  KINETIC  TEHP  (KELVIN)  IS  11:*/TK 

URXTE(6/*)  *THE  NCLECULAR  WEIGHT  (G/HOLE)  IS  12:*/XN 

ENDXF 

WRITE(E/*)  ‘THE  UPPER  STATE  VIBRATIONAL  LEVELS* 

WRITE(6^*)  *ARE  IN  EQUILIERIUP  (C>  OR  THEY  ARE* 

HRITE(6^*)  *EXCITEO  (1)  12:*/IfEXCXT 

URXTE<6^*)  *XNCLUDE  R-CENTRCIDS^YES (1 )»NO(C)  14:*/XFCENTR 

WRITECE^*)  ****UHICH  OF  THE  ABOVE  VARIABLES  DC  YOU*# 

1  *  WANT  TO  CHANGE?  TYPE  0  WHEN  YOU  ARE  FINISHED*# 

2  *  MAKING  CHANGES.  ENTER  NUMBERS  ONE  AT  A  TIME.* 

DC  ISC  I  =  1#U 

REAO((#*)  IVAR 
IF(IVAR.EC.C)  GO  TO  liC 
IF(IVAR.EG.I)  THEN 
ViRITE(6#*)  ‘PLOT  TITLE:* 

REAO(6#35)  NTITLE 
ENOIF 

IF(IVAR.EC.2)  THEN 
WR1TE(6#«)  *OISTHI8UTION  (0  OR  1):* 

REAO(E#*)  XFCIS 
XF(XFDIS.EQ.C)  THEN 

WRXTE(6#*>  *REMEM3ER  TO  MODIFY  VARIABLES  11  AND  12* 

ENOIF 

ENDXF 

XF(IVAR.EC.3)  THEN 

bRITE(6#*)  *ENVE10PE  HALF  WIDTH  (ANGSTROMS):* 

READ(6#*)  ENVHI 
cNDIF 

IFdVAR.EG.A)  THEN 
bRXTE(6#«)  *EFANCHIN6  IS;* 

REA0(6#*)  IBRANCH 
ENDXF 

XFdVAR.EG.S)  THEN 

WRITE(6#*)  ’MIN  WAVELENGTH  (ANGSTROMS):* 

READ(G#*)  WMIN 
ENDXF 

IFdVAR.EQ.G)  THEN 

kRITE(6#«)  *MAX  WAVELENGTH  (ANGSTROMS):' 

READ(6#*)  WMA> 

ENDXF 

XFdVAR.EO.?)  THEN 

WRXTE(6#<»)  *RESOLUTION  (ANGSTROMS):* 

REA0(6#*)  RRES 
ENDXF 

XFdVAR.EO. 8)  THEN 

WRXTE(6#*)  *Vie  TEMPERATURE  (KELVIN):* 

REAO(G#*)  XTV 
ENOIF 

IF(IVAR.EC.9)  THEN 

WRXTE(6#*)  *RCT  TEMPERATURE  (KELVIN):* 
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REA0(6^*)  XTR 
END  IF 

IF(1V*R.EG.10)  THE^ 
feRITE(6^'*)  IS:* 

READ(6/*)  6AFFA2 
ENOIF 

IFdVAfl.EC.II)  THEN 

fcRlTE<6^*>  *CAS  KINETIC  TEMPERATURE  (KELVIN):* 

REAO(6/*)  TK 
ENOIF 

IF(IVAR.EC.12)  THEN 

URITE(6^*)  *FCLECULAR  WEIGHT  ( GR A M S /MC L £  )  : * 

R£A0(6/’*)  XM 
ENOIF 

if<ivar.ec.13)  then 

kRlT£(6/>*)  *LF»ERSTAT£  VIERATIONAL  EXCITATION* 
fcRIT£(6#0  *ECUrLIBRrLW  tC)/  NON- E  OU  U  I BR  1 UM  (1):* 
flEAD(6^<-)  IFEXCIT 
END  IF 

IF(IVAR.EC.IA)  THEN 

WRITE(6^*)  *1NCLUCE  R-CENTRCICS?  TESd)  OR  NO(C):* 

REACCE/*)  IFCENTR 
ENOIF 

150  CONTINUE 
160  CONTINUE 

IFdFEXCIT.ES.I)  then 
CALL  ASSIGN(2/ *PCPLTN.WT*  ) 

REA0(2r*)  NVIE 

C  NVIE  IS  THE  NUNEER  CF  VIBRATIONAL  LEVELS  IN  THE  UPPER  STATE  YOU  HAVE 
C  WEIGHTS  FOR.  ALSO  THERE  IS  ONE  SET  OF  WEIGHTS  FOR  EACH  TRANSITION. 
00  6C  N=1/NF 

REA0<2/*)  <WGHT(N/ri)/II*1.NVlB> 

WHITE  (3^*)  (WeHT(N/II)/Il*1,>NVIB) 

60  CONTINUE 

CALL  CLCSE(2) 

ENOIF 

IFdFCENTR.EQ.I)  THEN 
CALL  AS3IGN(2. 'RCENTRCIO.DAT') 

00  85  N=1.NF 

C  N=1  IMPLIES  BI/2  —  Xl/2  ELECTRONIC  TRANSITION 
C  N=2  IMPLIES  83/2  --  X3/2  ELECTRONIC  TRANSITION 
C  THERE  SHOULD  BE  AS  PANT  R-CENTROIOS  AS  FRANK-CCNOON  FACTORS 
00  9C  II=1/NF1 

R£AD<2/*)  <CNTR<N^ri^IP)^IF=1^NF2) 

WRITE  (3/-*)  (CNTR  (N/  II^IP)/^IP=1  rNF2) 

90  CONTINUE 
35  CONTINUE 

CALL  CL0SE(2) 

ENCIF 

WR1TE(6^*)  *INCLUDE  NON-BOLTZMANN  CONTRIBUTION  TO  ROTATIONAL* 
WRITE(6/*)  'LEVELS  AS  SPECIFIED  8ELCU  ?  YES  (1)/  NO  (0).* 
READ(6/*)  IFPEPT 
IF<I FPERT.EQ.1)  THEN 

WRITE(6#*)  *WHAT  IS  THE  FSEUOC  RCTATICNAL  TEPFERATURE  (K)?* 
READ<6^*)  TRP 

WRITE(6x*)  *WHAT  IS  THE  KELAT  HEIGHT  CF  2N0  EXP  FOR  POT  CIS?* 
READ(6^*)  VAR 

WRITE(6/<‘)  'AT  WHICH  ROT.  LINE  OCES  CCNT.  TAKE  EFFECT* 

REA0(6^*)  VA«1 

VAR2SVAR1-1 
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ENOIF 

WRITECE^*)  •PUBLICHTION  QUALITY  (GRAPHS  BLRCK)?  YES  (Dr  NC  (0).* 
RER0(6r*>  IFPUE 

WRITECEr*)  *00  YOU  WANT  TO  PLOT  THE  THECRETICAL  ANB  EXPERIPENTAL * 
bRITECEr*)  'SPECTRA  ON  THE  SAME  GRAPH?  YES  (D  r  NC  (0)* 

READCEr*)  IFDATA 
IFdFOATA.EQ.D  THEN 

WRITECEr*)  'SPECIFY  THE  THE  AVERAGE  NUPEER  OF  EACKGROUKC  CCLNTS* 
WRITE(Er*)  'PER  PIXEL  TC  SUBTRACT  FROM  YOUR  DATA.  DC  THIS  BY' 
WRITElEr*)  'FIRST  SPECIFYING  THE  INITIAL  NIPBER  OF  COUNTS  TO' 
URITE(Er*)  'SUBTRACT  AT  YOUR  MINIMUM  WAVELENGTH' 

REA0(6r«)  EKGRNC 

WRITE(Er*)  'NOfc  SPECIFY  THE  SLCPE  CF  THE  LINE  CF  THE* 

URlTECEr*)  'AVEFAGE  INCREASE  IN  COUN TS / ANST R OM ' 

READCEr*)  SLOPE 

HRITE(Er*>  'YOL'f  SPECTRAL  DATA  IS  ASSUMED  TC  BE  IN  THE  ' 
WRITECEr*)  'FILES  NAMED  IN  :  FILE.CMA.  NOTE  THAT  IT  IS' 
WRITE<dr*)  'ALSC  ASSUMED  THAT  THE  FIRST  FILE  LISTED  IN  FILE.OMA' 
WRITE(Er*)  'CORRESPONDS  TO  THE  LOWEST  WAVELENGTH  AND  SC  ON  ' 
WRITECEr*)  'IT  IS  ALSO  ASSUMED  THERE  ARE  NC  GAPS  IN  WAVELENGTH* 
WRITE(Er*)  'BETWEEN  CONSECUTIVE  OMA  FILES.* 

ENOIF 

IF(I  FDATA.NE.DTHEN 

Wfl:TE<Er*)  'CALCULATING  THE  SYNTHETIC  SPECTRA...' 

WRITECEr*)  ' CRUNCH... CRUNCH. ..CRUNCH...' 

ENDIF 

ENOIF 

TVsXTV*.695 
Tfl*XTR*.E95 
IFCIFOIS.EC.D  THEN 
V»SQRT<2.49AEa*TK/XM)*1£8 
ENDIF 
RETURN 
END 

SUBROUTINE  INTENSE 
PARAMETERCNR0T»125) 

C0MHCN/FACT/F(3r30r 30)rYC9,30r2)r6C9r3C)r 
1  8(9r3C)r06(9r3C) 

COMMON/P  /TRW ( 6 r 3 Or 3C )r Pi < 6r KROT ) r P2 CEr NRCT) r 

1  P1INTC6rNROT)rP2INTCErNROT) 

COMMCN/C  /01CErNR0T)rQ2(ErNRCT)r 
1  Q1INT(6rNRCT)rC2INT(6rNR0T) 

COMMON/ R  /Rl(ErNR0T)rR2CErNRCT)rR1INT(6rRR0T)r 
1  R2lNT(ErNfiCT)rFINTC6r1Q00C)rFINC1CC00)rFIMXC1CC00) 
COMMCN/TEMP/TVrTR 

COMMON/ INPlT/NTITLErlFDISrlERANCHrWMINrWMAXrRRESr 

1  XTVrXTRrIFOONEr If DEF^IMINrlMAXr lENVHr ENVHirVr 

2  XMrTXr6AMMA2 

COMMON/ 2NOEX/NFrNF1rNF2rNTrN0rN01 

C0MMCN/W6HT/W6HTC3Cr30)rNVIGrIFEXCIT 

C0MM0N/CNTR/CNTR(3r3Cr3C)rIFCEKTR 

COMMON/EXPD AT/I FO AT Ar WAVE (20 rlC30)r SPEC (20r103C)rWAV(2CC00)r 

1  SP(2CCCO)rWW(2COOC)rSPW(2CCCO)rWOATC2C)rSENC2C)rlCCUNTrIIMAXr 

2  INMAXrBKGRNDrSLOPE 
COMMON/SPIN'/NC 
CHARACTER*4C  NTITLE 


C 

C 

C  CALCULATE  ENERGIES  (G)  AND  KCLECULAR  CONSTANTS  (E)  FCR 

^5 


C  VIBR•TIO^AL  LEVELS 
C 

NM/IX=N«X(NF1  /NF2) 

DO  10 

00  20  I*1/NrAX 
G(f»/I)  =  0. 

DO  30  L^UNDI 
C 

IEP=L-1 

XI*FLOXT<I) 

6(P^ I)*G(W#1>*Y (M^L^I )*(<XI-1.)*.5)**(I£P) 

B(P^I)  =  B(Pl^I)^YCMi-L#2)*  ((X1-1.>4.5)**(IEP) 

C  WPITE<3#*)  ‘G/P^rs'/GCP^D^W^I 

30  CONTINLE 

WR1TE<3#*)  *G/P/1=*#G<P#I)/K^I 
C 

C  CALCULATE  THE  ENERGY  DIFFERENCE  COG)  BETWEEN  VIBRATIONAL 
C  LEVEL  I  AND  THE  L0'.4EST  VIBRATIONAL  LEVEL  WITHIN  THE  ELECTRONIC  STAT 
C 

DG(H^I)=G(P^I)-G(M^1) 

HRITE<3^*)  ‘DG/P^I^'rOGCN^D^M^I 
C 

20  CONTINLE 
10  CONTINLE 
C 

C  CALCULATE  THE  VIERATICNAL  LEVEL  TRANSITION  (TR) 

C  FRECUENCIES  EETWEEN  THE  UPPER  AND  LCWER  ELECTRONIC 
C  LEVELS  --  HINDS  CASE  A. 

C 

NNsI 

DO  40  Nal^NF 
DO  50  I  =  nNF1 
DO  60  IP=1/NF2 

TRW(N^I^IP)»G<NN*1/I)-G(NN#IP) 

WRITE  <3  MRANS*/TRW<N^I»IP>^N<-I^IP 
C  NOTE:  EP-BPP<C  FCP  RED  DEGRADED  BANOS 
C  SP-BPP>C  FCF  BLUE  DEGRADED  BANDS 

WRITE<3/*)  *BP-BPP»*^B<NNt1,I>-B(NN^lP> 

60  CONTINUE 

50  CONTINUE 
NN=NN>2 

40  CONTINUE 
C 

C  DETERPINE  THE  PIN  AND  MAX  EXTENT  OF  THE  PESH  IN  ANGSTOMS 

C  IF  NOT  SET  IN  THE  INPUT 

C 

IF(IFDEF.EG.O)  THEN 

WM1N-1E20 

WMAX>C. 

DO  41  N»1^NF 
DO  51  l=nNF1 
DO  61  IP*1^NF2 
XW»1.E8/TRW(N^1/IP) 

WMAXsMAXCXW/WMAX ) 

WHIN-MIN(XW^WMIN) 

61  CONTINUE 

51  CONTINUE 

41  CONTINUE 
RRES>(WP.AX*WMIN)/2C00. 
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ENOIf 

WRITE(3#*)  •RES**^PRES/  'WMIM**  *  WR#X»%  WHIJi 

IF(ENVRI.EC.C.C)  E^VH1«RRES«40. 

IERVH*2NT(ENVH1/RRES) 

IMINsIM<CWriN^.1*RRES)/RRES)-IENVH 
IH/^X*I^T<(UR/^X4.1*flRES)/RRES>♦IENVH 
WRITE (3^*)  •ENVM  =  */ENVH1^*I«IN»*/IPIIN/*I«AX  =  *  ^IHAX 
C 

C  READ  IN  THE  EXPERIMENTAL  WAVELENGTHS  AND  PHCTON  COUNTS  IF 

C  THE  SYNTHETIC  SPECTRUM  IS  TC  BE  COMPARED  WITH  EXPERIMENTAL 

C  RESULTS- 
C 

IFCIFDATA.EQ.I  )  THEN 
CALL  LCAOFIL 

WRITE(6#*)  'SPLICING  AND  CALIBRATING  THE  EXP  SPECTRAL  FILES...* 
CALL  SPLICE 

CALL  hEIGHTCwMIN/VMAX) 

WRITECE/*)  'CALCULATING  THE  SYNTHETIC  SPECTRA...* 

WRITECE/*)  'CRUNCH... CRUNCH... CRUNCH...* 

ENDIF 

C 

C  FOR  EACH  VIERATICNAL  TRANSITION  CALCULATE  THE  MLLTITtDE  CF 
C  TRANSITIONS  CUE  TC  ROTATIONAL  SPLITTING  CF  THE  VIBRATIONAL 
C  LEVELS  AND  CALCULATE  THE  RELATED  INTENSITIES. 

C 

WMlNE  =  taM.IN-2*ENVHl 
WMAXEsfcMAX*2*ENVH1 
NN^  2 

DO  75  N*1/NF 

C  THE  EOUILI0RIUM  POPULATION  CF  NS1/2  AND  NSJ/2  LEVELS 
C  ARE  NOT  THE  SAME.  THE  RELATIVE  NUMEER  CF  MCLECLLES  IN 

C  EACH  LEVEL  IS  EASILY  DETERMINED  USING  THE  ELECTRCNIC 

C  DISTRIBUTION  FUNCTICN.  AT  3CO  DEGREES  KELVIN  IT  SHOWS 
C  THAT  IF  THE  POPULATION  OF  NS1/2  IS  GIVEN  BY  1/  THEN  THE 

C  POPULATION  CF  NS3/2  IS  GIVEN  BY  .6  .  THEREFCRE^  THE 

C  N'S  IN  THE  HERZEERG  INTENSITY  FCRMULA  DIFFER  AS  MENTIONED 
C  ABOVE  BETWEEN  THE  1/2  AND  3/2  LEVELS. 

IFCN.EC.DTHEN 
FELEC  =  1  . 

ELSE 

FELEC=.6 

ENOIF 

C  DETERMINE  THE  NORMALIZATION  FACTOR  FCR  BCLTZMAN 
C  EQUILIBRIUM  DISTRIBUTION  OF  POPULATIONS  IN  THE  UPPER 
C  STATE 

Z  =  C. 

DO  70  JJ*1^NF2-1 
Z*Z+EXF (-DG (NN/ JJ)/TV) 

70  CONTINUE 

00  80  I-1«NF1 
DO  90  1P*1/NF2 
WIL>1.E8/TRk(N^I^IF) 

IFCWIL.LT.WMINE.OR.WIL.GT.WMAXE)  GO  TO  90 

c 

c  THIS  IS  THE  NCRFALIZEO  POPULATION  FCR  EACH  VIBRATIONAL  LEVEL 

C  IN  THE  UPPER  ELECTRONIC  STATE 

C 

IFdFEXCIT.EQ.C)  THEN 

EXPV*EXP(-0G(NN/I>/TV)/2*FELEC 

WGHT(N#I)«EXPV 
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ENCI  F 


C 

C  IF  THERE  IS  A  K  CN-E  flUIL  IB  RIUF*  CISTRIBUTION  CF 

C  VIERATIONAL  LEVELS  THEY  ARE  SET  HERE 

C 

IF(IFEXCIT.E0.1 .ANC.NVie.GE.I)  E XPVs U6HT CN^ I ) 

C 

C  IF  ONE  REQUIRES  THE  INCLUSION  CF  THE  R-CENTROICS  IN  THE 

C  COPPUTATION  OF  THE  INTENSITIES  THEIR  CCNTRIBUTICN  IS 

C  TAKEN  INTO  ACCCINT  NEXT 

C 

IFdFCENTR.EC.I)  THEN 

RE  =  -T.  +  1.76  711«CNTR<N/I/>IP)-1.C375  5*CNTR<N/I^IF)**2 
1  ♦2.C28*CNTR<N/1/IP)**3 

EXPV=EXPV*RE**2 
ENDI  F 
C 

C  THE  LAMEDA  TYPE  COUELING  OF  EACH  MUL T IPLET (TWO  CF  THEM  NC»1/2) 

C  IS  TAKEN  INTC  ACCOUNT  IN  THE  NC  LOOP.  THE  EEFF  CF  THE  FUNCTIONS 

C  FP  AND  FPP  CHANGE  ACCORDINGLY. 

C 

DO  72  NC=1/2 

IF(I6RANCH.EQ.1 )  CALL  D CU9L E Tf N, NN , I ^ I P # E XP V ) 

WRITE (3^*)  *F/I^IP/EXPV**^F<N#I^IP>/I/IF/EXPV 
C 

C  REDISTRIBUTE  INTENSITIES  WITH  LCRENT2IAN  “DUE  TC  HOMOGENEOUS  SLIT 
C  RELATED  EROADENING 

C 

C  IF(IFDIS.EC.O)  THEN 

CALL  LORN72N<N/I/IP) 

72  CONTINUE 
90  CONTINUE 
€0  continue 
NN=NN+2 
75  CONTINUE 
RETURN 
END 

SUE  ROUT INE  0 OUE L ET ( N / NN / I ^ IP ^ EXP V ) 

PAflAMETER<NROT*125) 

C0MMCN/FACT/F<3»304  30)#VC9#30^2)#-G(9<3C)^ 

1  9<9^3C)^D6<9^3C) 

COMMCN/P  /TRW<6#30/30)^P1 <G#NROT) ^P2 <C ^NROT) # 

1  P1INT(6#NRCT)^P2INT<6^NR0T) 

COPMON/G  /01(6zNR0T)#C2Ce^NRCT)^ 

1  Q1INT(6/NRCT)^C2INT(64fNROT) 

COMMON/ R  /R1 <E/NR0T)#R2<6,NRCT)/R1 INT(6.NR0T)# 

1  R2INT<«^NRCT)/FINT<6/1C0CC)^FIN(1CCCC>^FIKX(1CCC0) 
COMMON/TEMP/TV/TR 

COMMCN/INPUT/NTlTLE^IFOIS,ieRANCH#kMlN/ViMAX#RRES^ 

1  XTV/XTR#IF00NE«IP0EF/IPIN#IMAX/IENVH#ENVH1#V# 

2  XM^TK^6AM,MA2 

COMMON /IN0EX/NF/NF1/NF2/NT^ND#N01 

DATA  0E#0EP^AP#APP/1 .2E-6^1.3E-«#149.7/285.e/ 

C  YY  AND  YE  ARE  THE  LOWER  AND  UPPER  STATE  SPIN  ORBIT  SPLITTING 
C  WAVE  NUMBERS. 

DATA  P#PE^YY/YE/1./1./221.5#90./ 

C  DATA  P/PE#YY^YE/1.^1.^100CC^1CC0C/ 

CHARACTER*4C  NTITLE 
IF(N.EC.I)  THEN 
2*. 5 
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ZE».5 

ELSE 

Z-1.5 

ZEsl.5 

ENDIF 

DO  100  J«1/NR0T 
IF(N.EC.I)  THEK 

C  N  E8  1  P'E/^^S  UE'RE  C0I^6  BP  11 /2  — XPIl /2  TRAKSITIONS 
C  PI-ERAhCH 

R*FL0/iT(J>-.5 

RE*R-1 

PI  (N^J)*TPW<N^1^IP)<'FP(R-UNN/I^P/YE/0E)-FPP<B^NN/  IP^F^  YY/OEP) 
Ffl=P1<N#J) 

C  NOTE:  THE  CCNPILER  bAS  HAVING  A  PRCBLEP  RECOGNIZING  Pi  IN  THE 
C  XI  FUNCTION  CALL  UNTIL  I  SET  FRsP1(iy/>J)  ANC  PASSED  FR  IN  THE 
C  XI  FUNCTION  CALL/  AS  YCU  SEE  BELOW.  BEFORE  THE  CHANGE/  AT  RUN 
C  TIPE/  FR  WITHIN  THE  XI  FUNCTION  CALL  WAS  UNDEFINED  »  ZERO. 

CALL  AINT<SJ/F/R/2/YY/RE/PE/ZE/YE) 
P1INT(N/J)»XI(N/NN/I/IP/SJ/FR/RE/EXPV/YY) 

C  WRITE<6/*)  ‘SJ  OF  P1»»/SJ 

C  WRITE<6/*)  N/ J/Pl {N/J)/P1INT(N/J) 

C  QI'ERANCH 
R£*R 

Q1 <N/J)sTRW<N/l/ IP )+FP<R/NN/I/P/ YE/DE )-FPP(R/NN/IP/P/YY/OEP) 
FRsOI (N/J) 

CALL  AINT(SJ/F/R/Z/YY/RE/PE/ZE/YE> 
01INT(N/J)«XI(N/NN/I/IP/S4/FR/RE/EXPV/YY) 

C  R1-BRANCH 

RE»R41. 

R1  < N/ J  )»TR  W<N/1/ IP  )->-FP{R^1/NN/ I/P/ YE/DE  )-FPP(R/NN/ IP/P /YY/OEP) 
FR«R1 (N/J) 

CALL  AINT(SJ/P/R/Z/YY/RE/PE/ZE/YE) 
R1INT<N/J)*XI<N/NN/I/IP/SJ/FR/RE/EXPV/YY) 

ELSE 

C  N  EC  2  MEANS  WE*RE  DOING  8P 1 3/2--XP 1 3/2  TRANSITIONS 
C  P2-eRANCH 

R-FLCATCJI+.S 
RE*R-1 . 

P2<N/J)*TRW(N/I/IF)+FP<R-1/NN/IP/P/YE/DE)-FPP(F/NN/IP/F/YY/DEP) 

FRaP2<N/J) 

CALL  AINT(SJ/P/R/Z/YY/RE/PE/ZE/YE) 
P2INT(N/J)«X1(N/NN/I/IP/SJ/FR/RE/EXPV/YY) 

C  WRITE(6/«)  ’SJ  of  P2**/SJ 

C  WRITE(6/*)  N/J/P2(N/J)/FR/P2INT(N/J} 

C  02-eRANCH 
RE»R 

Q2(N/J)»TRH(N/I/IF>*FP(R/KK/IP/P/YE/0E)-FPP<R/NN/IP/P/YY/CEP) 

FR>Q2(N/J) 

CALL  AINT<SJ/P/R/Z/YY/RE/PE/ZE/YE) 

Q2INT (N/J )«XI(N/NN/ I/IP/ SJ/FR/RE/EXPV/YY) 

C  R2>ERANCN 

RE«R*1. 

R2 ( N/ J >sTR W (N/ 1/ IP ) ♦FP(R^1/NN/1P/F/ YE/C E)-FPP(P/NN/ IP/P /YY/DEF) 
FR»R2(N/J) 

CALL  AINT(SJ/P/R/Z/YY/RE/PE/2E/YE) 

R2INT(N/J)«XI(N/NN/I/ IP/ SJ/FR/RE/EXPV/YY) 

ENDIF 

100  CONTINUE 
RETURN 
END 
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RE»L  FLNCTICN  F F < R/ NN # I E #DE ) 

PABAKElEREKROTsIZS) 

CCPPCN/FAC7/FC3,30-r30>^Y<9^3C^2)/G<9/3C)#- 
1  9<9/3C>/D€(9/3C) 

COPMCN/F  /TR^<6^30/3C)#P1<«,^B0T)^P2(<^^NRCT)/ 

1  P1INT(6/NBCT)/F2IM(6/NR0T) 

COPMON/SPIK/NC 
IF(P.LI.O)  BETIBN 
IF(NC.E0.1>7HEN 

BEFF*B<NN/T)*<1.*B(NN^I)/(F*A8S(Y£>)) 

ELSE 

BEFF  =  B<NN^1)*(1 . - B ( NN# I ) / < F * ABS ( Y E )  )  ) 

ENCI  F 

FP=BEFF*R*(fl*1)-DH*R**2*(R^1)**2 

RETURN 

END 

REAL  FL'KCTICN  F  F  P  (B IP  ^P  ^  Y  Y#  D  EP  ) 

PARAPETER(hROT=125) 

COPMCN/ FACT/FC3/3 3/ 30)^7(9/30, 2)^G(9/3C)/ 

1  3(9/3C)/D6  (9/30 

COPPCN/F  /TRfc<6/30/3Q)/Pl <6/NR0T) /P2 (E/NPCT) / 

1  P1INT(6/NBCT)/F2IKT(6/AR0T) 

COPKCN/SPIK/NC 
IF(B,L7.C)  BETLBN 
IF(NC.EC.1)TH£N 

B£FF=e(KN-1/IP)*(1.*B(NN-1/IP>/(P*AES(YY)>> 

ELSE 

eEFF=E(KN-1 /IP)*(1.-9(NK-1/IP)/(P*AES(YY))) 

ENCI  F 

FPF  =  B£FF*R*(R  +  1)-OE*R**2*<P-»1)**2 
C  IF(N.E0.2}  THEA 

C  WBITE<6/*)  *N/KN/!/B<NA/r>/e<AN-1/IP>/F/PF* 

C  UBlTE(6/«)  N/KK/I/6<NN/I)/E(NN*1/IP)/B/PF 

C  EKOIF 

RETURN 
END 

REAL  FLNCTICN  > I ( N/ NN/ I / 1 P/ S J / FR / R, E XP V / Y Y  ) 
C0PP0N/FACT/F(3/30/ 30)/Y(9/30/2)/G(9/3C>/ 

1  3(9/3C)/0G(9/3C) 

C0PHCN/7EPF/TV/7R 

C0PMCN/PER7/IFF£RT/TRP/IFPUe/VAR/VAB1,VAR2 
DATA  P/DE  /1/T,2E-</ 

IF(B.GE.O)  THEN 

X1=B<NK/I)/TR*SJ*FR**3*F<N/I/IP>* 

1  EXF(-3(NN/1)*R* (R+1 )/TB)*EXPV 

ELSE 
XI=0 
ENDIF 

C  DUE  TO  PERHAPS  A  SLCULY  RELAXING  NASCENT  DIET,  C«  PERHAPS 
C  DUE  TO  CURVE  CBCSSINGSC?)  AT  R  VALUES  GPEATER  THAN  APPRCXIPATEIY 
C  50  (SEE  JEFFERIES)  IN  V*  =  C/1  WE  SEE  A00I71CNAL  PCFULATIONS 
C  IN  THE  HIGHER  BCTATIONAL  LEVELS.  THIS  EFFECT  IS  APFBOXIPATED  BY: 
IFdFPERT.EC.DTHEN 
IF(I.LE.2.AND.B.G£.VAR1 )  THEN 
XI»XI  +  VAR*e(NN/l)/TRP*Sj*FI>**3*F(N/I/IP)* 

1  EXF(-B(NN/1)*(B-VAS2)*(B-VAR2+1 )/TRF)*EXPV 

ENDIF 
ENDIF 

C  IF(N.EQ.2)  then 

C  HBITE(6/*)  •N/NN/I/9(NN/I)/SJ/FR/F(N/I/IP)/EXPV/R' 
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C  URITE(6/*}  N^KR/I/B(NK«I)/SJ^FR/F(K#I/IP)/EXPV#R 

C  ERDIF 

RETURN 
ENC 
C 

SUBROUTINE  AINT ( SJ « P/ R« 2/ V# R E« PE ^ ZE # Y E ) 

C  THIS  SUBROUTINE  CAtCULATES  THE  ROTATIONAL  LINE  STRERGTHS  FOR 
C  DOUBLET  TRASNSITIONS. 

C  DEFINITION  OF  SYMBOLS  IN  AINT  SUBROUTINE: 

C  INPUT  CONSTANTS  <LCWER  STATE  IS  A  SINGLE  LETTER:  EXCITED  STAT 
C  HAS  AN  ADDED  E): 

C  RxRE  -  ROTATIONAL  QUANTUM  NUMBER/  J  (POSITIVE  HALF  INTEGER). 
C  SELECTION  RULE:  OJ=0/1/-1 

C  P/PE  -  angular  momentum  CF  ELECTRONIC  STATE/  LAMCA  (POSITIVE 
C  INTEGER  OR  ZERO).  SELECTION  RULE:  OLAMDA=C/1 /-I 

C  Z/ZE  -  TOTAL  AN6LLAR  MOMENTUM  OF  THE  ELECTRONS/  CME6A 
C  (HALF  INTEGER  *  LAMOA  ♦OR  *  .5) 

C  Y/YE  -  DIMENSIONLESS  SPIN  ORBIT  COUPLING  C ON S T AN T / L A MD A. 

C  PURE  CASE  A/  LAMDAsICCOO.  PURE  CASE  B/  LAMDA<E-5. 

C  OUTPUT: 

C  SJ  -  ROTATICNAL  LINE  STRENGTH  FACTOR 
C 

SJ=O.C 

SL=P-PE 

IF<AeS<SL)-1-01)1/1/35 

1  SRsR-RE 
IF(ABS(SR)-1. 01)2/2/35 

2  IF<ABS(Y)-«999. 05/5/3 

3  IF<ABS(YE)-9999. 0)13/13/4 

4  IF<AeS(2-ZE)-. 01)13/13/35 

5  IF<ABS(Y)-. 00002)4/6/13 

6  IF<ABS(YE)-. 00002)7/7/13 

7  IF<A8S<R-2)-. 01)10/10/8 

8  IF<ABS<RE-2E)-. 01)10/10/9 

5  IF<A6S<SR-SL+Z-ZE)-1.01)10/1C/35 

10  IFCP-, 01)11/11/13 

11  IFCPE-. 01)12/12/13 

12  IF<A8S(SR'AZ-ZE)--01)35/35/13 

13  IF(R-Z4. 01)35/14/14 

14  IF<RE-ZE*. 01)35/15/15 

15  AJ»RE+1.0 
SsP-.S 
SE«PE-.5 

16  IF<ABS<SL)-. 01)23/17/17 

17  IF<A6S<SR)--C1)18/21/21 

18  Q».5*SCRT((AJ^SL*SE)*CRE-SL«SE) 

1  •(SE^AJ)/(RE*AJ)) 

GO  TO  27 

19  Q».5*SCRT((AJ«1.«SL*SE)*(AJ«SL*SE)AAJ) 

GO  TO  27 

20  Q«SE*SCRT((AJ4RE)/(RE*AJ)) 

60  TO  27 

21  IF(SR)22/19/19 

22  0*.5*SCRT((RE-SL*S)*<AJ-SL*S)/RE) 

60  TO  27 

23  IF(ABS(SR)*. 01)20/24/24 

24  IF(SR)  26/26/25 

25  0«S0RT((AJ*AJ*SE*SE)/AJ) 

GO  TO  27 

26  G»SORT((RE*RE>S«S)/RE) 
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27 


28 

29 

30 

31 

32 

33 


34 

35 


1 

2 

3 

4 

5 


6 

7 


8 

9 

1C 


C 


SI»1  -0 
SK«1 .0 
SIE«1,C 
SKE«1 .C 

I F CABS <Z-S)-. 01 >30/30/28 
SI»-1.C 

IF<Z-S“».01)29/29/3C 

SK»-1-C 

IF(ABS<ZE-SE)-.C1)33/33/31 

SIE=-1.C 

IF<ZE-SE+.C1)32/32/33 
SKE=-1 .C 
S  =  S+  1 
SE=SE+1 

SJsSJ+7<Sl/SK/P/R/\)*9*1<SlE/SKE/PE/RE/1E) 

IF<S-P-1-49> 18/34/24 

SJ=SJ*SJ 

RETURN 

END 

REAL  FLNCTICN  T < SI/ SK/P/R/ Y) 

IF<P-.C1>8/8/1 
1F(R+.49-P)9/9/2 
AJ  =  2. *P  +  1  . 

IF(AeS(Y)-.C00C2)8/8/3 

SL=Y/AeS<Y) 

IF<ABS<Y)  *9  999. 04/4/7 

gasy*p*p/aj-r-.5 

68s.5*SCRT<AJ**2+Y* (Y-4.)*F*P) 

BK»1 ./SCRT(2.*6E*AJ) 

T*e»:*<SCRT<<6e-SI*EA)*(R*.5-Sl.*P>)^Sl 

1^S0fiT<(G9+SI*GA)*<F+.5^SL*P)>>*SK 

GO  TO  1C 

TsSI/SflRT(2,0)*SK 
GO  TC  10 
GA=SL*P*P/A J 
G8=- 5*F 
GO  TO  5 

T*SI*SA*SCRT(<R+.5-SI*P)/AJ) 

GO  TO  1C 
T*.5*SK«<1 .♦SI) 


RETURN 

END 

SUBROUTINE  LCRNT ZN(N/ 1/ IP > 

PAPANETER<AROT*125) 

C0PM0N/FACT/F<3/ 30/ 30)/Y<9/3C/2)/G<9/30)/ 

1  9(9/3C)/DG  (9/30 

COPMCN/P  /TRV<6/30/.3C)/P1(6/NROT)/P2(6/NROT)/ 

1  PlINT<6/NRCT)/F2lAT(d/AR0T) 

COPNCN/G  /Q1<8 /NR0T)/fl2(8/NRCT)/ 

1  01INTC6/NROT)/C2INT(6/NROT) 

COPNON/R  /R1C8/NR0T)/R2«/NRCT)/R1INT(6/NR0T>/ 

1  r2INT<8/NRCT)/FINTC0/1C00C>/F1N(1CCCC>/FIA)K1CCC0) 


1 

2 


COPPCN/TEMP/TV/TR 

COPNON/INPUT/NTlTLE/IFOIS/IBRASCH/VPIN/Vi-AX/RRES/ 

XTV/XTR/IFOCNE/IFCEF/IPlN/IPAX/IEKVM/EhVHl/V/ 


XP/TK/6APPA2 

COPFICN/INOEX/NF/NF1/NF2/NT/AD/N01  _ 

C0PNCN/PEA*{/XPEAX(8/30/30)/YPEAK  (6/30/30  /FCHCICOOC) 


CHARACTER*4C  ATITLE 
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NORFALIZEO  LCREATZIAN— SLIGHT  R E D I  ST R IB 17 1  OK  OF 
INTENSITIES 


C 
C 

c 

c  6/IMHA2s3.16ET3«3,16E13 

IOIF=IFAX-IHIN 
IXN=3 

WIL=1-E8/TRH<N/I/IP) 

IF<WIL.LT.WHIN4ENVH)  THEN 

IL*INT((HIL'».1*FRES)/RRES)-IHIN*1 

IU*IL*IENVH 

XN^WIL 

ELSE 

IL=INT((U1L'».1'‘FRES)/R<<ES)-IHIK-IENVH^1 

IU=IL+2*IEKVH 

XNaWIL-ISNVH*RPES 

ENDIF 

WRITE  (3/*)  •  IL  =  '/IL/*IL**#IU 

FMAX»C 

FMAX1=C 

FnAXCLC*0 

FMH  =  C 

RR  =  C. 

00  110  I1=U/-IL 

IF(I1.LT.1.CR.I1.GT.I0IF)6C  TO  110 

RRsRR+RRES 

XI1=XN+RR 

C  XI1I=3.£ie/XI1 

SUH»G 

IFCN.EC.I)  THEN 
00  120  J*1/NRCT 

0EN0Hs<XI1-1E8/Fl (K  ^ J )) •*2*6 AMHA2 

SC*6Af1FA2/DeNCF 

SUF=SUF  +  SC*PHKT  <N^  J  ) 

120  COKTINLE 

00  130  J=1/NR0T 

0EKCHs<XI1-1  .E€/R1<N^J))**2-»GArf»A2 

3C=GAHFA2/0ENCF 

SUF*SUF  +  SC*PHKT<N/J) 

130  CONTINIE 

DO  140  J«1/NRC7 

0EK0r=<XI1-1.E?/Q1  (Nz-J)  )**2^6AHHA2 

SCs6AnPA2/CEN0H 
SUH>SUP‘»SC«C1  IKT  (N#  J) 

140  COKTINLE 
ELSE 

DO  150  J=1/KPCT 

DEK0«»(XI1-1  .E8/P2<N/J)  )**2'»GAHf'  A2 

SC*GAMPA2/CEN0H 

SUHsSUH-»SC*P2lNT  <N/J  ) 

150  COKTINLE 

00  160  J*1/NRCT 

0EK0H*<XI1-1-Ee/R2<N/^J))**2'*GAPHA2 

SC=6AHPA2/OENOP 

SUraSUP^SC«R2lK7 (N/ J) 

160  COKTINLE 

DO  170  J«1/KRC7 

0EKCM=<XI1-1  .E8/e2(N^J)  )**2  +  6AF'HA2 

SC«GAnPA2/CENCP 

SUP«SUP^SC*C2IKT (N/J) 

170  COKTINLE 
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ENOIF 

FIN(I1)=FIN(I1)*SUM 
FCH(I1 )=SUF 

C  WRITE(6/*)  •0EN  =  ’/0ENCf1^SC^6#f»MA2»Q2UT(N/J)#SUM 

C  FIND  THE  FAX  T  PCSITICN  OF  "CLOSE*’  PEAKS 
FPAX  =  PAX<FIN(I1)/F(1AX) 

FPAX1=PAX<FCH<I1)#FPAX1) 

IF ( FMAX.GT.FMAXOLC . AND. (FPAX1 /FKAX ).GT..2)  THEN 
FPAXCLD=FPAX 

YPEAK  (N/1/IP)  =  PAX<YPEAK  <K/.IyIP>^FFAX) 

END  IF 

11C  COMINLE 
C 

C  END  CF  IKTEGPATICN 

C 

C 

C  FIND  THE  PEAK  OF  EACH  SPECTPUP  LINE 
C 

FMAX=0 
FMAX1=C 
FMAXCLD=C 
RR  =  0. 

IL=INT((WU>.1*Pfi£S)/RRES)-IF:iK-IXK^1 

IU=IL+2*IXK 

XN=WIL-IXN*RR6S 

DO  5  ri=IL^IU 

IFdl.LT.UCR.n.GT.  IDIr)6C  TO  5 

RRsRR^PRES 

XI1=XN4RR 

FFAX=FAX(FIN<I1)/FMAX) 

FPAXlsPAX(FCH(I1)/FHAX1) 

IF<FMAX.GT.FHAXCLC,ANO.<FPAX1 /FHAX) .GT..2)  THEN 
XPEAK  CN/I/IP)  =  XI1 

YPEAK  (N,1^IP)  =PAX  <YPEAK<N^I/IP)/FPAX) 

FNAXCLD=FPAX 

ENOIF 

5  CONTINIE 
RETURN 
END 

SUBROUTINE  G AUS S N <N # I ^ I P) 

PAPAP,ETER(NR0T  =  125) 

COPHCN/FACT/fC3/30/30)^Y(9/30r2)/G<9#3C)^ 

1  B<9/3C)^DG<9^3C) 

COPKCN/P  /TRN<6/30^30>/P1 <6^NROT)rP2(C^NPOT)^ 

1  P1INT(6^NRCT)/F2INT(6^NR0T) 

COYMON/C  /Q1 (d/NROT),C2CE^NRCT)^ 

1  Q1INT<6/NRCT>/C2INT<6^NR0T) 

CCPPCN/P  /R1  <d/-NROT)^R2(6/.NRCT)^R1  INT  (6/NROT)  # 

1  R2INT<6/NRCT),PINT<6,1COOO^FIN(1CCOC)^FINX<1CCCO) 
COPKCN/TEMP/TV/TR 

COPPCN/INPLT/NTITLE/IFOIS^ISRANCH/kPIN/WHAX/RRES^ 

1  XTV/XTR/IFDONE/IFOEF^IPlN/iMAXi-IENVH/ENVHl^V# 

2  XP^TK#6AMPA2 

COPPCN/INOEX/NF/NFl/NF2^NT/NO/NOl 

COPPCN/FEAK/XPEAK<d/-30^30)#YPEAK(6/3Ci’3C)/FCH(1CCOC) 
CHARACTER*4C  NTITLE 
C 

C  NCRPALI2EO  GAUSSIAN 
C 

RETURN 
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END 

SUBROUTINE  OUTFIT 
PARA«ETER<KR0T*125) 

C0FHCN/FACT/F{3#30/ 30)^T<9#30/2)^G<9^30)< 

1  8<9/3C>^0G<9#3C) 

COFNCN/P  /TRVi(«/3  0#30)i-P1  <6^PROT)/P2<e/NRCT)# 

1  P1INT<6/NR0T>/F2INT(6/KR0T) 

COKMCN/G  /Q1 (E/NR0T)«02(6/NR0T># 

1  Q1INT<6^NROT)/C2IAT(6^NROT) 

COPMCN/R  /R1 (E/NR0T),P2C6^NRCT)#P1INT«6#KR0T)/ 

1  R2lNT(6/NFCT>/FINT<6/1C00C)/fIK(1CCCC)/FIN>(1CC00) 
COPMCN/TENP/TV/TR 

COPMCN/IN'PLT/NTITLE^IFOIS^IEPANCH^WPIN/IiNAX/PRES^ 

1  XTV^XTfi/IFOCNE/IFDEF^IPlN/lPAX^IEAVM/EAVHWV/ 

2  XP/TK^6ARPA2 

COPKCN/INOEX/NF/NF1 ^NF2^NT^KO^KDl 

COPM0N/PEAK/XPEAiC(i/30#3C)/VPEAX(6^20/>3C)xFCH(1C00C) 
COPWCN/WGHT/UGHT  (3C^30)/NVI6r  IFEXCn 
C OP rCN/C NT fl/C NT F  (3/30/30  ^IFCEKTR 

COPMCN/EXPCAT/IFOAT A/WAVE<2C/1C30)/SPEC(20/103C>/WAV(2CC00)/ 

1  SF (20CCO)/WU(2COCC)/SPN<2CCOO)/WDAT(2C>/SEN(2C)/ICCUNT/IIPAX, 

2  IKMAX/EK6PN0/SL0PE 

C0PPCN/FERT/IFFERT/T9P/IFPLE/VAF/VAP1/VAR2 
OIPENSICN  F  INYdOCCO)  /  YKAXXC 30/30) 

CHARACTERoAC  NTITLE 


CALL  ASSIGN (7/*CUTFUT. SPEC*  ) 


WRITE  INPUT  FILE  FOR  CONTINUATION  RUNS 

CALL  ASSIGN<2/* INPUT. OAT* ) 

WRITE<2/*)  WPIN/WHAX/RRES/XTV/XTR 

URITE<2/35)  NTITLE 

F0F«AT(A3C) 

WRITE  (2/*)  IF0IS/IERANCP/ENVH1/TN^XP/GA»'"IA2 
WRITE<2/*)  IFEXCIT/IFCENTR 
CALL  CLCSE(2) 

N0RPALI2E  INTENSITY 

YMAXsO- 
YniN*1.E20 
LHNalNAX-If IN 
LMNN=LPN-2*IENVK 
DC  10  lal/LMNN 
00  1C  I*IENVH/LPN-IENVH 
YnAXaMAX(FIN(I)/YnAX) 

YMIN=NrN<FINCl)/YNIN) 

COKTINLE 

YM*Y«AX 

II«0 

WNENaUPIN-ENVNi 
DC  20  lal/LNNN 
00  20  1*IENVH/IPN-IENVH 
II*IIt1 

FINX(Il)aWPEN^FLOAT (I-1)*RRES 
FINV<II)»FIN<I)yVMAX 
UR1TE<7/*)  FINXdD/FINVCiI) 
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20  CONTINtE 
C 

C  N0RFALI2E  PEFK  IMEKSITIES 
C 

00  25 
DO  24 

00  23  IP=1/NF2 

YPE4K<K/I/IF)syFEAKN^r#IP)/VM/‘X 

23  CONTINUE 

24  CONTINUE 

25  CONTINUE 

YMIN  =  YPIN/YPA)( 

YM4X=1 . 

C 

C  PLOT  THE  INTENSITY  SPECTRUM  AS  A  FUNCTION  OF  WAVELENGTH  IN 
C  ANGSTROMS  USING  CISSPLA. 

C 

TTV=TV/.695 

TTfl=TR/.695 

C  XLWAV=WMIN-ENVM 

C  XLkAV1=WMAX>ENVHl 

XLkAV=:kPIN 
XLWAV1=WMAX 
CALL  PAGE(11^8.5) 

CALL  AREA2C(1Q/£.5) 

NCHAR=INOEX(NTITLE/*  *) 

CALL  HEAOIN  <XREF (NTITL£)/NCHAR^2.^1 ) 

CALL  XNAMEC'UAVELENGTH  < ANG STRCM S  )  S * /I CC ) 

CALL  YNAMECNORPALIZEO  INT£NSITYS»/1CC) 

CALL  GRAFCXLWAV/’SCALE'  /XLkAV1^YMIN^*SCALE%YMAX) 

CALL  MESSAG<*RCT  TEMPaS  *#-1CC/7.35#-6.25) 

CALL  REALNC{TTR/-4/*AB*JT*#*ACUT») 

CALL  DASH 

CALL  CLRVE<FINX/FINY/LMNN/C) 

CALL  RESETC'OASN*) 

WRITE<3/*)* IF0ATA^IFEXCIT=«^IF0ATA/IFEXCIT 

IFCIFEXCIT.EQ.OTHEN 

NVI0=T2 

CALL  MESSAGC’Vie  TEMP=$V1CC/7. 35^6.0) 

CALL  REALNC<TTV#-4/»ABUT*^*AeUT*) 

IF(IFCATA.EQ.C)  THEN 

CALL  MESSAGC’V  1/2-1/2PCF  3/2-3/2P0P! * CO/7 . 35^ 5 .75 ) 
YPOS»5.5 

CO  2CC  JJ>1^NVZ9 
XP0S=7.35 

CALL  INTNC<JJ-1/XPCS/YPOS) 

XPOS=XPCS*,3125 

CALL  REALNOCkGHTCI^JJ )/-2^XPCS^YPCS) 

XP0S=XP0S-*1.3 

CALL  REALN0ChCNT(2#JJ )^-2#XPCS^YPCS) 

YP0S=YP0S-.2 
200  CONTINUE 
ELSE 

CALL  PESSAGC*--  SYNTHETIC  SPECTRUPS* #100^7. 35#5. 75) 

CALL  PESSAGCV  1/2-1/2PCP  3/2-3  A  2PCPS  *  / 1  C0#7 . 35/ 5. 25  ) 
YPCSaS.O 

00  210  JJa1/NVIB 
XPCS*7.35 

CALL  INTNC<JJ-T/XPOS/YPOS) 

XP0SaXP0S4.3125 
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CALL  AEALN0(heHT(1^JJ)#-2#XPCS^YPCS) 

XP0SSXP0S41.2 

CALL  PEALN0(hCHT(2/JJ }#-2#XPCS#YP0S) 

YPCS»YP0S-.2 
210  CONTINUE 

IFdFPUe.EQ.C)  THEN 
CALL  SETCLRCAEO') 

ENDIF 

CALL  NESSAGC*  -  EXPEP.  SP £ CTROPS* / 1 00/ 7. 2 « / 5. 5 ) 

C  CALL  PARKERCIS) 

C  CALL  SCLPIC(.7) 

C  CALL  CURVE(bL/SPW/INPAX/INPAX/1CC) 

CALL  CURVE(Uk/SPU/INPAX/C} 

WRITE<7/*)  'EXP  UAVELENSTH  /  EXP  FHCTCN  COUNTS* 

00  27  I1>1/INPAX 
b9ITE(7/«)  Wb<ri)/SPWCII) 

27  CONTINUE 

EKOIF 
ENCIF 

IF(IFEXCIT.EO.I)  THEN 
IF  (IFCATA.EQ.C)  THEN 

CALL  PESSAGCV  1/2-1/2PCP  3/2-3/2PCPI * /1Q0/7. 35/6.0) 

YPOS*5.75 

CO  3CC  JJ«1/NVIB 

XP0S*7.35 

CALL  INTACCJJ’I/XPOS/YPOS) 

XPOSSXPOS+.3125 

CALL  PEALNOCWCHTd  /JJ)/-2/XPCS/YPCS> 

XP0S«XP0S41.2 

CALL  PEALNC<b€HT<2/JJ)/-2/XPCS/YPCS) 

YP0SSYP0S-.2 
200  CONTINUE 
ELSE 

CALL  PESSAGC*--  SYNTHETIC  SPEC  TRIPS */ 100/7 . 35 / 6 .C ) 

CALL  PESSAG(*V  1/2-1/2PCF  3/2-3/2P0FJ * / ICO/7 . 35/ 5- 5 ) 

YPOS»5.25 

00  31C  JJsl/NVIB 

XPCS»7.35 

CALL  INTAC<JJ-1/XP0S/ VPOS) 

XPOS=XPOS>*.3125 

CALL  REALN0{b6HT(1/J4)/-2/XPCS/YPCS) 

XP0S*XP0S-H,3 

CALL  REALKC{VCHT<2/JJ)/-2/XPCS/YPCS) 

YP0S»YP0S-.2 
310  CONTINUE 

IFCIFPUB.EO.OTHEN 
CALL  SETCLPCPEO*) 

ENOIF 

CALL  PESSAGC*  -  EXPEP-  SPE CTRUPS * / 1 OC/7. 3 5 /5. 75 ) 

C  CALL  PARKER(IS) 

C  CALL  SCLPIC(.7) 

CALL  CURVE(yb/SPb/lNPAX/C) 

bRITE(7/*)  'EXP  bAVELENGTH  /  EXP  PHOTCN  CCUNTS* 

00  127  II>1/INPAX 
bRITE(7/*)  WbCID/SPWCII) 

127  CONTINUE 

ENOIF 
ENOIF 

CALL  NElGHTC.Ce) 

CALL  AN6LE(90.) 
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XN=1 .7*.007353*<«M4X-WHIN) 

YN=.025 

CALL  SETCLR (*BLLE*) 

C  FIND  UHEffE  LINE  TRAKSITICN  NLMBERS  WILL  CVERLAP  AND  ELLIPNATE  THE 
C  CVEBLAPPING  NUHEERS 
N  =  1 
NN  =  2 

DO  4C  I=1^NF1 
DO  50  IP=1/NF2 

IF(AES(XPEAK(N^1/IP)-XPEAK(KN/1/IP>).LT.XN>THEN 
WRITE!  I/ •)  •ABS(XP1/2-XP3/2)yXN/IxIF='rASS(XFEAK(N^l#IP)- 
1  >PEAK(NN^I/IP))^XN^I/IP 

YPEAK(NK/1»IF)=YPEAX(NN#I#IF)*1C 
WRITE(3i-*)  •YPEAX(2/’I/IF>=*.rYPEAK(NN/I/’IP) 

ENDIF 
50  COMINLE 
40  CONTINUE 

C  WRITE  THE  NC K-OV E 9LA oP IN G  NLP8EPS  CN  THE  PLCT 
N=1 

DC  60  I=l^Nf1 
DC  70  IP  =  UNF2 

IF(XLWAV.GT.XPEAK(K,I,IP).Cfi.XLWAVl.LT.XPE*K(N/I/IP)) 

1  GC  TC  70 

1F(YPEAX(N^I/-IF).LT..001)  GC  TC  7C 
IC= 1-1 
IF0=IF-1 

XPESS=XPEAX(N/I/IF) 

YH£SS  =  1.t-YN 

CALL  FLNESSC  <  %1/XWeSS/YPESS) 

CALL  RLINT<IO/ ’ABLT'/'ABUT*) 

CALL  fiLPESS<*/*/1/  '  A8UT*^*Ae*JT*  ) 
call  BLINT(IPC/*AEJT'^*ABtT’) 

CALL  RLMESS <•) 1/21*/ IOC/ • AEUT»#  ’AEUT* ) 

70  CONTINUE 
60  COMINLE 
N*2 

DC  00  I=1/NF1 
DC  90  IP=1/NF2 

IF(XLUAV.GT.XPEAK(N,I,IP).CP.XLWAV1.LT.XPEAK(N/I/IP)) 

1  GO  TC  90 

IF(YPEAK(N/I/IP).LT..001)  6C  TC  90 
IF<YPEAtC(N/I/IP).LT.8)  THEN 
ICaI-1 
IPC=IP-1 

XHESS3XPEAK(N/I/IF) 

YPESS*1.*YA 

CALL  RLMESSC*  (  ' /1/XF'ESS/Y»'ESS) 

CALL  RLINT(IC/*ABLT*/*A0UT*) 

CALL  fiL«ESS<*/*/1/*ABUT»/*AeOT*  ) 

CALL  BLINT(IPC/*AEUT'/*ABLT*) 

CALL  RLFIESS(*)3/2S*/1CC/*AEUT*/*AELT') 

ENOIF 

90  CONTINUE 
00  CONTINUE 

CALL  RESETCHEIGHT*) 

CALL  RESETCAKGLE*) 

CALL  EN0PL(-1) 

C 

RETURN 
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subroutin;  PLTOEV  (MAHDEV) 

(OEFIKE  SUBROUTINE  HERE  !!!) 


WRITTEN  EY: 


FRANK  J.  SEILER  RESEARCH  LA6CRATCRY 
UNITED  STATES  AIR  FORCE  ACAOERY 
COLORADO  SPRINGS/  COLCRAOO  8C84C 


INPUT  PARAMETERS: 

NAMOEV  -  OUTPUT  DEVICE  CHOSEN  BY  USER 


CHARACTER*4S  NAMDEV 


*  LOCAL  VARIABLES: 

*  I0UF 

*  lERR  -  ERROR  NUMBER 


INTEGER  IBUF/I/II/IERR 

DIMENSION  IBUF(16) 


CHARACTER  EELL 

BELL  s  CHARC?) 

*  TERMINAL  DEVICES 

IF  (NAMDEV  -EC.  'VT  240')  GOTO  21C0 
IF  (NAMDEV  .EC.  'TEKTRONIX  41U7'>  6CTC  22CC 

IF  (NAMDEV  .EC.  'TEKTRONIX  4C10')  GOTC  23C0 

IF  (NAMDEV  .EC.  'VT  ICO  RETRO  GCTO  240C 

*  HARDCOPY  DEVICES 

IF  (NAMDEV  .EC.  '  PRINTROMCS' )  GCTO  31CO 
IF  (NAMDEV  .EC.  'HP  7550')  GCTO  320C 

IF  (NAMDEV  .EC.  'HP  7475')  GCTO  3300 

IF  (NAMDEV  .EC.  'ANY  ASCII  PRINTER')  GOTC  54CC 
IF  (NAMDEV  .EC.  'POST  PRCCESSOR')  GOTC  41C0 

*  ERROR/  QUIT 
GOTO  910C 

*  TERMINAL  DEVICES 


*  VT  24C 

2100  CONTINUE 

*  SET  CCNFIGURATICN  TO  STANDARD  I/O 

IBUF(1)«C 

CALL  I0MGR(ieLF/-102) 

CALL  VT24G 
GOTO  990CC 

*  TEKTRONIX  41C7 

2200  CONTINUE 
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*  SET  CCNFICL'RATION  TO  STANDARD  I/O 

I0UF<1)=C 

CALL  I0F'GR<I8LF/-102) 

CALL  TK41  (A1C7) 

*  PLOT  ORIENTATION  ’ 

* 

* 

C  CALL  HWICTC'CCF'IC') 

6CT0  99CCC 

*  TEKTRONIX  4C1C 

2300  CONTINUE 

CALL  P401C 
GOTO  9?CCC 
VT  1C2  UITH  GPAPMCS  RET^CFIT 
240C  CONTINUE 

*  SET  CCNf IGURATICN  TO  STANDARD  I/C 

*  I0UF(1)=C 

CALL  I0f<GR(IELF/-1C2) 

CALL  PDEVT 
SOTO  59CCC 

* 

2500  CONTINUE 

GOTO  990CC 


2600  CONTINUE 

GOTO  990C0 

* 

2700  CONTINUE 

GOTO  990CC 
*  HARDCOPY  DEVICES 
«  PRINTRCNICS 
310C  CONTINUE 

CALL  PPNTNX 
GOTO  99CC0 


AlTO  :  ALTO  ACCORDING  TO  PLT  SIZE 
CCPIC:  X  AXIS  HORIZONTAL 
PCVIE:  Y  AXIS  HORIZONTAL 
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*  HP  7550 


3200  CONTINUE 

*  SET  CCNFIEURATICN  TO  SPOOL 

IBUF(1)-S 

CALL  I0KCR<IELF/-102) 

*  SET  FILE  OPEN  MODE 

*  3  -  NO  OVERWRITE^  CREATES  NEW  VERSION  NtPBER 

*  C  -  APPEND  TO  FILE 
IBUF<1)=3 

CALL  lO^CRClEt F^-1C4) 

SET  PLOT  FILE  NAPE 

00  3220  11*1^16 
IBUF(II)>C 
CONTINUE 

PRINT  *<////)* 

CONTINUE 
PRINT 

••  ENTER  PLOT  FILE  NAPE  OR  0  TO  QUIT:  BELL 

READ(5#  '(UAA)*/  ERR*93CC)  (IBUF(II)/  11*1/  U) 

IF  (leUF(l)  .EQ.  ) 

CALL  I0P.Gfi<ieLF/-1C3) 

CALL  HP755Q  <1) 

PLOT  ORIENTATION  -  ALTO  :  AL'TO  ACCORDING  TO  PLT  SIZE 

CCPZC:  X  AXIS  HORIZONTAL 
PCVIE:  Y  AXIS  HORIZONTAL 

CALL  HHRCTC’PCVIE*) 

GOTO  99CCC 

*  HP  7475 

3300  CONTINUE 

*  SET  CONFIGURATION  TO  SPOOL 

IBUF<1)»5 

CALL  ICHGR<IELF/-102> 

*  SET  FILE  OPEN  PCCE 

*  3  -  NO  CVERWRITc/  CREATES  NEW  VERSION  NLPBER 

*  0  -  APPEND  TO  FILE 
IBUF(1)*3 

CALL  I0HGR<IELF/-1C4) 


*  SET  PLOT  FILE  NAPE 

DC  3320  11*1/16 
ieUF(lI)«G 
3320  CONTINUE 

PRINT  •<////)• 

3340  CONTINUE 

PRINT  •<••!••/ 

♦  **  ENTER  PLOT  FILE  NAME;  "/A)*/  BELL 


3220 

3240 


* 

* 

C 


REA0(5/  CRR  =  94CC)  (I9UF(II)/  11  =  1^  14) 

CALL  ICMCR <IEL F/-103) 

CALL  HP7475  (1) 

PLOT  ORIENTATICN  -  ALTO  :  ALTO  ACCCROIMC  TO  PLT  SIZE 

COPIC:  X  AXIS  hOBIZOMAL 
POVIE;  Y  AXIS  HORIZOKTAL 

CALL  HWRCTCAOVIE*  ) 

GOTO  99CCC 


•  ANY  ASCII  PRINTER 


* 

* 

* 

C 


3400  CONTINUE 

CALL  FPRTPL 
GOTO  990CC 
*  POST  PROCESSOR 
4100  CONTINUE 

CALL  COWPRS 


GOTO  99CCC 

•  ERROR  HANDLING 


*  DEVICE  SELECTION  ERRCR  - 

9100  PRINT  '(**0**/ 

*  «•  ***  ERROR  IN  PLOT  DEVICE  SELECTION  ***•*, 

*  A^  A/  A)*^  8ELL^  SELL#  SELL 
STOP 


*  HP 

9300 

4 

♦ 


7550  FILE  NAME  ERROR  - 

PRINT  *<*'0**# 

21X#  ••  ***  ERROR  IN  FILE  NAME  ***••#  A#  /# 

21X#  •*  PRESS  RETURN**#  A#  A)»#  EELL#  BELL#  EELL 

READ<5#  *(I1)*#  I0STAT*IERP)  I 
GOTO  324C 


*  HP  7475  FILE  NAME  ERROR  - 

9400  PRINT  *('*0'*# 

♦  21X#  **  *•*  ERROR  IN  FILE  NAME  ***•*#  A#  /# 

♦  21X#  *•  PRESS  RETURN**#  A#  A)»#  EELL#  BELL#  EELL 

READ(5#  *(I1)*#  IOSTAT=IERR)  I 

GOTO  334C 


99000  CONTINUE 
RETURN 
END 
C 

c «*«****«*«««««*««***««******«**«•**«***«***««**«««*«««*« «**««**«•••««**< 

SUBRCUTINE  LCADFIL 
C 

C  This  LOADS  THE  OPA  FILES  THAT  INCLUDE  WAVELENGTH  INFCRHATION 
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C  IT  ALSC  places  THE  WAVELENGTH  AND  CO P a ES PCND IN G  INTENSITY  DATA 

C  IN  FILES  TO  EE  MANIPULATED  IN  SUBROUTINE  SPLICE.  IT  IS  ASSUMED 

C  THAT  THE  CMA  FILES  TC  BE  PLOTTED  ARE  PRESENT  IN  THE  LOCAL  FILE 
C  SPACE  AND  LISTED  IN  FILE  FILE. CMA. 

C 

COMMON/ EXPC AT/I F DAT A, WAVE <20/1 C30) /SPEC (20/1 030  rWAVCZCCOO)/ 

1  SP<20CC0)/WU(2C00C)/SPfc<2CCC0)/HDAT  (2C)/SEN<2C)/ICCUNT/IIMAX^ 

2  INMAX/BKGRND/SLOPE 
COMMCN/COF/ ACOF  (20/50) 

CHARACTER<*6  SHOTNO/  SFC 
CHARACTER*^  END 
CHARACTEP-1Q  SNOCAT/FILES 
CHARACTER*2C  CDUMMY1..C0JMY1 
C 

c 

DUMMY=0 

ICOUNT=0 

FILES=' FILE. CMA' 

CALL  ASSIGN(A/FILES) 

liRITE<6/*)  'LCADING  FILES  (NAMED  IN  FILE-CMA):' 

DO  ICC  Ia1/2C 

fiEA0(4/1C1/END=9)  SHCTNO 

:C0UN7*IC0UNT41 

WfiIT£<6>.*)  I/'  FILE=  '/SHCTNC 

INUM=I 

101  FORMAT  <A6) 

END** .CMA  * 

SNOOAT*SHOTNO//END 
CALL  ASSIGN<2/SNCDAT) 

R£A0<2/1C2)  CCUMY1 

102  F0RMAT<A19) 

DC  1C  11*1/3 
fi£AD(2/1C3)  CCUMY2 

103  FORMAT(AII) 

10  CONTINUE 

DO  2CC  J=1/1024 
PEAO<2/*/ENO=2C1)  SPEC(I/J) 

200  CONTINUE 

201  CONTINUE 
PEAD(2/1C2)  CCUMYI 
DC  2C  11*1/3 
READ(2/1C3)  CCUMY2 

20  CONTINUE 

DC  21C  J*1/«0 

R£AD<2/*/ENO=211)  ACCF(I/J) 

210  CONTINUE 

211  CONTINUE 
CALL  CL0SE<2) 

C  WRnE<3/*)  'EXP  WAVE  LEN  /  EXP  PHOTON  CCLNT' 

DC  22C  J*1/1C2« 

WAVE(I/J)sACCF(I/36)+AC0F(I/37)*(FL0AT(J)-1.) 

C  WRITE<3/*)  WAVE<I/J)/SPEC<I/4> 

220  CONTINUE 

100  CONTINUE 

9  CONTINUE 

CALL  CLOSECA) 

RETURN 

END 

C 

SUERCUTINE  SPLICE 


43 


C  THIS  SUfiRCUTIHE  SPLICES  TOGETHER  CWA  FILES.  IT  IS  RSSL'MEC 
C  THPT  THE  CRTA  FILES  OVERLAP  EACH  OTHER  AND  THAT  THEY  ARE 
C  NAHEO  IN  FILE  FILE.CMA,  IT  IS  ALSO  ASSUMEC  THAT  THE  FIRST  FILE 
C  NARED  CORPESPCNDS  TO  THE  SHORTEST  WAVELENGTH  AND  SC  ON. 

C 

COHMCN/E)(PCAT/IFDA7A/WAVE(2C^1C30)/SPEC(20>103C)^WAV(2CCOO)^ 

1  SP(20CCQ)^UW(2CCCC}^SPk(2CC00)/UDAT(2C)/SEN(2C)/>ICCUNT/IIPAX, 

2  INMAX/BKCRND/SLOPE 
COHNCN/COF/ACCF (20^50) 

DC  20  I  =  1xICCL’M 
IFCI.EQ.D  THEN 
IC  =  1 

X=ACCFCl/36) 

ELSE 

lC=ICCLO-INT(<WMAXOLD-ACCF(I/36)>/DXCLO) 

X  =  X-(V*f«AXCLD-ACOF  (1x36)) 

ENOIF 

WRITE <3x*)  •ACOF(Ix36)=*xACOF{Ix36)x'ACOF(Ix77)='xACCF(Ix37) 
WRITE(3x*)  '  IC/XxWH!AXCLOxICOLO/DXCLD  =  *xICx  X/ VPAXCLOx  ICOLDx 
1  OXCLO 

OC  30  11=1x1024 
fcAV(IC>=X 

SP{IC)=HAX(S?(IC)xSPEC(IxII)) 

X=X*ACOF<!x37) 

IC=IC*1 
30  CCNTINL’E 

WH AXCLO  =  X-ACC  F  (Ix  37) 

1C0LC=IC-1 

CXCLC=ACCF<Ix37) 

20  CONTINUE 

IIPAXsICCLO 

END 

SUBROLTINE  WEIGHTCWPINxWMAX) 

C 

C  THIS  ROUTINE  WEIGHTS  THE  SPECTRAL  COUNTS  WITH  THE  INVERSE  OF 

C  THE  CPA  SENSITIVITY.  IT  ALSC  NOPKALIZES  THE  EXPERIMENTAL  SPECTRA 

C  AND  LOADS  THE  FINAL  PLOTTING  ARRAYS. 

C 

CCPMCN/EXPCAT/IF0ATAxWAVE(2C^1C3O)xSPEC(2Ox103C)xWAV(2CC00)x 

1  SPC20CCC)xWW(2CC2C)xSPW(2CC00)^W0AT(2C}xSEN(2C)xICCUNTxIIPAXx 

2  INPAXxBKGPNDxSLOPE 

DATA  WCAT  !  2 5C . x 26 0- x 270 2 80.x 29C . ^ 3CC - x 3 2 0. x 3 50. x 37C . x4CC . / 

1  4  5O^.x5C0.x5  5O-x6CO.x6  5C-x7CC.x75  0.x800.x90C.x1C00/ 

C  DATA  SEN  /  2 . 6 9/ 3 . 1 6x 3 . 32 x 4. 26 x 5 . 2 2x4 . 2x 5 . 5 3 x 5 . 55 x 5 . C4 x 6 . 1 x 

C  1  8. 24x11. 73x17. 81x29. 7x49. 7x57. 03x117. x237.x513.7x 

C  1  7E7.97/ 

DATA  SEN  /  .529x.576x.6x.634x.672x.755x.?56x1.C6x. 972x1. CCx 

1  1.35x1.94x2.99x4.93x8.0x9.5x19.34x37.93x1038x1669/ 

C 

UUPIN=WPIN/10. 

WyPAX=WPAX/10. 

11  =  0 

SPHHAX=1E-2C 

WRITE(3x*)  •lIPAX»*xIIPAX 
00  100  I=1xIIHAX 

c  yRiTE(3x*)  'wwpiNxyyNAXxyAv=*xwyniKxywPAXxkAV(i) 

IF(WAV(I).LT.yWPIN.0R.yAV(I).6T.yyPAX)  GOTO  110 
J*0 

II«II+1 

IF(yAV(I).LT.W0AT(1))  THEN 
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«T»SEK(1) 

ELSE  IF(WAV(I>.LT.IiOAT(2>)  THEK 
J*2 

WT*SEI»<J-1)*(WAV(l)-W0AT<J-1))«(SEh<J)-SEK<J-1))/ 
1  <H0AT(J)-WCATCJ-1)) 

ELSE  ZF(UAV(lX.LT.hDAT(3))  THEN 
J*3 

WT*SEA<J-1)^<LAV<I)-W0AT(J-1) )*<SEh<J)-SEA<J-1))/ 
1  (WDAT(J)-W0AT<J-1 )) 

ELSE  IFCWAV<I>.LT.UDAT(4))  THEK 
Ja4 

UTaSEK(J-1)+(UAV<I)-W0AT<J-1)>*<SEK(J)-SEK(J-1))/ 
1  <W0ATCJ)-W0AT<J-1)) 

ELSE  IF(WAV(IX.LT,W0AT(5))  THEK 
JaS 

WT«SEK  (J-1)*(ViAV<I)-M0ATCj-l)  )*(SEN(J)-SEK(J-1  ))/ 
1  (WCAT(J)-W0AT<J-1  )) 

ELSE  IF(UAV{I).LT.kDAT(6))  THEK 
J  =  6 

WT*SEK<J-1)*(iiAV<I)-W0AT<J-1))a(SEK<J)-SEK(J-1>)/ 
1  <yCAT<J)-WCAT(J-1)) 

ELSE  IF  (HAV <I).LT.fcOAT<7) )  THEK 
J*7 

WTaSEK (J-1 )^(KAV<I)-W0AT<J-1) )*(SEK (J>-SEK (J-1 >)/ 
1  <kCAT<J)-WDAT<J-1 )) 

ELSE  IF(WAV(I).LT.VOAT(e))  THEK 
Js8 

WTaSEK<J-1)+(KAV<I)-U0AT<J-1))*(SEK<J)-SEK(J-1))/ 
1  <UCATCJ)-kOAT(J-1)) 

ELSE  IF(HAV<I).LT.kOAT<9))  THEK 
J*9 

WT*SEK<J-1)^<UAV<1)-W0AT<J-1)>*(SEK(J>-SEK(J“1))/ 
1  (W0AT<4)-H0AT<J-1)) 

ELSE  IF  (UAV(I)«LT.kOAT(10))  THEK 
JalO 

WTsSEK<J-1)  +  <»iAV<I)-UDAT(J-1))*CSEK(J)-SEN(J-1 ))/ 
1  <WCAT(J)-HOAT<J-1)) 

ELSE  IF<HAV(I).LT,V0AT<11))  THEK 
Jail 

HT«SEK<J-1)^(KAV<l)-HOATCJ-1))a(SEK(J)-SEK(J-1))/ 
1  <H0AT<J)-y0AT<4-1)) 

ELSE  XF(WAV(I).LT.VDAT(12))  THEN 
Ja12 

HT*SEK<J-1)>(KAV<I>-HOAT<J-1))a(SEK<J)-SEN<J-1))/ 
1  <WOAT(J)-«OAT<4-1)) 

ELSE  IF(WAV(I).LT.K0AT(13))  THEK 
Ja13 

WTsSEK(J-1>+(hAV<I)-M0AT<4-1)>aCSEK(J>-SEK(J-1))/ 
1  («CAT<J)-HCAT<J-n) 

ELSE  IF(UAV(I).LT.KOAT(14))  THEK 
J*14 

HT*SEK<4-1)*(kAV(I>-HOAT(J-1>)a(SEK<J)-SEK<J-1))/ 
1  <HOAT(J)-«DAT<J-1)) 

ELSE  IF(yAV(I).LT.K0AT(15))  THEK 
Ja15 

HT«SEK<J-1)^CfcAV<I)-H0A7(J-l))*<SEK<J)-SEK<J-1))/ 
1  <«CAT<J)-HDAT(J-1>) 

ELSE  lF(WAV(X).LT.K0AT(16n  THEK 
Ja16 

WTaSEK{J-1>+CWAV(I>-HOAT<J-1>)*<SEK<J)-SEK<J-1))/ 
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1  (UD/IT(J)-WD#T(J-1)) 

ELSE  rF<UAV<I).LT.hOAT(17))  THEN 
J  =  17 

WT=SEN  (J-1  )*(liAV(I)“H0ATCi-1))*CSEN  (  J  )  -  S  EN  (  J -1  )  )  / 

1  (WCAT(J)-UC»T(J-1)) 

ELSE  IF(WAV(I).LT.ViOAT(18))  THEN 
J*18 

WT=SEA(J-1)*(fcAV(I )-WDAT(J-1))* (SEN (J)-SENt J-1  ))/ 

1  (WCAT(J)-WD«T(J“1)) 

ELSE  IF(WAV(I).LT.N0AT<19))  THEN 
J  =  1  5 

WT=SEN(J-1>^(fcAV(I )-WOAT(J-l) )* (SEN (J)-SEN(J-I))/ 
i  (WCAT(J)-liCAT(J-l)) 

ELSE  IF(WAV(I).LT.W0AT(20))  THEN 
J  =  20 

WT=SEN  (J-1  )  +  (l.AV(I  )-H0AT(J-1))*  (SEN  (J)-SES(J-1))/ 

1  (WCAT(J)-WCAT(J-I)) 

ELSE  IF (WAV (I) .GT.UOAT(20))  THEN 
WT=SEN(20) 

ENOl  F 

C  SHIFT  EXP.  SPECTRA  TO  LEFT  EY  6  ANGSTRCPS  OLE  EXP.  CAL.  ERRCR 
UW(II)=UAV(I)*1C.-6 

SPW(II)=SP(I)*h7-(6XGRNC^SLCP£*(WU(II)-WMIN)) 

C  SPW (II)=SF (I)*fcT 

C  ELLIPINATE  NH  PEAX  FFOP  DATA  N CRP A L I Z A T  I CN--C R  OTFER  PEAKS  IF 
C  SC  DESIRED. 

IF(3345.LT.WW(II),AN0.3375.GT.VW(II) 

1  .CR.525C.LT.WW(H).AN0.53AC.GT.WW(II))  THEN 
C  DC  nothing 
ELSE 

SPWNAX=PAX<SPW(II)^SPUrAX) 

ENtlF 

110  CONTINUE 
100  CONTINUE 
INPAX=II 

WRITE(3^*)  ' INPAX=», INPAX 

C  NORMALIZE  THE  EXPERIPENTAL  SPECTRA  TC  THE  HIGHEST  NON* 

C  EXCLUDED  PEAK 

DO  2C0  I  =  UINPAX 
SPW(I)=SPW(I)/SFUMAX 

C  HRITE(3/*)  •Wfc/SPW=*/HW(I)/SPW(I) 

200  CONTINUE 

C  CHOP  OFF  UNWANTED  LIKES  EXTENDING  OUTSIDE  OF  THE  PLOTTING 
C  AREA 

DO  21C  I=1/INPAX 
IF(SPW<I).GT.1.C1)SPW(I)=1. 

210  CONTINUE 
RETURN 


2  7  1C 

.C028  .0193 
.0137  .0658 
.C357  .1121 
.C656  .1238 
.C952  .0551 
.1164  .0484 
.1249  .0119 
.C031  .0211 
.C153  .0707 
.0395  .1176 
.0717  .1253 
.1027  .0508 
.1236  .0413 
.1300  .0072 


1.5951  1.6164  1 
1.5811  1.6C15  1 
1.5676  1.5877  1. 
1.5544  1.5740  1. 
1.5417  1.5606  1. 
1.5293  1.5472  1. 
1.5172  1.5331  1. 
1.5886  1.61C2  1. 
1.5744  1.5955  1. 
1.5606  1.5812  1. 
1.5472  1.5672  1. 
1.5342  1.5536  1. 
1.5215  1.540C  1. 
1.5091  1.5252  1. 


.0637  .1315  .1908  .2065 
.1329  .1365  .0608  .0013 
.1178  .C2C3  .0072  .0781 
.0458  .0048  .0738  .0557 
.0015  .0537  .0562  3.E-5 
.0143  .0695  .0043  .0412 
.0496  .0343  .0136  .0553 
.0679  .1369  .1944  .2061 
.1373  .1342  .0546  .0003 
.1151  .0246  .0111  .0837 
.0392  .0084  .0784  .0493 
.0002  .C6C5  .0501  .0011 
.0208  .0679  .0014  .0484 
.0573  .0270  .0206  .0517 


Franck-Condon  Factor®  (File: 


.1727  .1142  .0605  .0259 
.0344  .1171  .1606  .1375 
.0896  .0212  .0076  .0817 
4.E-5  .0511  .0849  .0245 
.0554  .0523  3.E-5  .0513 
.0467  .0006  .0564  .0383 
.0005  .0466  .0297  .0076 
.1691  .1096  .0570  .0239 
.0407  .1235  .1612  .1333 
.0855  .0159  .0119  ,0902 
.0003  .0585  .0820  .0181 
.0614  .0456  .0005  .0590 
.C4CC  .0021  .0613  .0304 
9.E-5  .0522  .0219  .0140 


FRANK. CON) 


6383  1.6607  1.6837 
'6230  1.6444  1.6656 
6080  1.6273  1.6588 
5930  1.6229  1.6379 
5699  1.6030  1.6216 
5715  1.5882  1.6005 
5567  1.5736  1.5997 
6322  1.6548  1.6780 
6168  1.6384  1.6596 
6017  1.6209  1.6513 
5864  1.6145  1.6316 
5473  1.5963  1.6151 
5636  1.5814  1.5881 
5492  1.5665  1.5918 


1.7074  1.7320 
1.6716  1.7186 
1.6758  1.6972 
1.6575  1.5644 
1.7599  1,6679 
1.6318  1.6498 
1.6155  1.6072 
1.7018  1.7265 
1.6460  1.7126 
1.669?  1.6914 
1.6512  1.7186 
1.6579  1.6616 
1.6251  1.6432 
1.6088  1.7C29 


1  .7410 
1.7410 
1.7169 
1.7078 
1.6872 
1.7060 
1.66C1 
1.7410 
1.7357 
1 .7103 
1.7017 
1.6809 
1.6818 
1.6534 


1.7A1C 
1.7410 
1.741C 
1.7290 
1.5298 
1.6977 
1.6771 
1.741C 
1.7410 
1.7410 
1 .7221 
1.7485 
1-6913 
1.67CC 


1.7410 

1.74tC 

1.7<1C 

1.741C 

1.7398 

1.7164 

1.7139 

1-7410 

1.7410 

T.741C 

1.741C 

1.7336 

1.7097 

1.7048 


R-centrolds  (File:  RCENTROID.DAT) 
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t 


0.  .7<9602 

1218-7  -.0081 

-7.28  0. 

30294.9  .5962 

797.31  -.0048 
-3.72  0. 

221.5  .775156 

1218.70  -.0063 

-7.28  0- 
3C384.1  .6013 

798.78  -.0046 
-3.59  0. 


Dunham  Coefficients  (File:  DUNHAM. COF) 


7 

-313  .110  .ICC  .C6  -C424  .0000625  .0CC003 
-188  -066  .060  .C360  .0254  -CC0037506  -OC0002 


Vibrational  Population  Weights  (File:  POPLTN.WT) 


f  3003C 

3203C 

3403C 

3603C 

38030 

40030 

4203C 

4403C 

46030 

48030 

5C030 

52050 


List  of  Experiment.-.!  OMA  Filec  (File:  FILE.OMA) 
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1.C0CC0C 


45C.00C0 


40CO-000  420C.C0C 

S<B->X)  SPtCl.U,,  .T  3 

1.00C0C0 

1  1 


12CC.CC0 

C.OOCCOOCE+OC  0 


Input  File  Generated  by  Code  (File:  INPUT.DAT) 


.COOCCOOE^CO 
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WfiVELENGTH  C ANGSTROMS) 

FIGURE  2 
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